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We present the results of a survey of 45 young (<250 Myr), close (<50 pc) 
stars with the Simultaneous Differential Imager (SDI) implemented at the VLT 



survey, we observed 54 objects total, of which 45 were close, young stars, 2 were 
somewhat more distant (<150 pc), very young (<10 Myr) stars, 3 were stars 
with known radial velocity planets, and 4 were older, very nearby (<20 pc) solar 
analogues. Our SDI devices use a double WoUaston prism and a quad filter to 
take images simultaneously at 3 wavelengths surrounding the 1.62 /xm methane 
absorption bandhead found in the spectrum of cool brown dwarfs and gas giant 
planets. By performing a difference of adaptive optics corrected images in these 
filters, speckle noise from the primary star can be significantly attenuated, result- 
ing in photon (and flat-field) noise limited data. In our VLT data, we achieved H 
band contrasts > 10 mag (5cr) at a separation of 0.5" from the primary star on 
45% of our targets and H band contrasts of > 9 mag at a separation of 0.5" on 
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80% of our targets. With this degree of attenuation, we should be able to image 
(5cr detection) a 5 Mj„p planet 15 AU from a 70 Myr Kl star at 15 pc or a 5 
Mjup planet at 2 AU from a 12 Myr M star at 10 pc. Our 45 southern targets 
were observed with the VLT while 11 of our northern targets were observed with 
the MMT (2 objects were observed at both telescopes). We believe that our SDI 
images are the highest contrast astronomical images ever made from ground or 
space for methane rich companions <1" from their star. We detected no tenta- 
tive candidates with S/N > 2 a which behaved consistently like a real object. 
FoUowup observations were conducted on 8 <2cr candidates (with separations 
of 3 - 15.5 AU and masses of 2-10 Mj^p, had they been real) - none of which 
were detected at a second epoch. In the course of our survey, we also discovered 
5 new close stellar binary systems with measured separations of 0.14" to 0.26". 
For the best 20 of our survey stars, we attained 50% 5cr completeness for 6-10 
Mjup planets at semi-major axes of 20-40 AU. Thus, our completeness levels are 
sufficient to significantly test theoretical planet distributions. From our survey 
null result, we can rule out (at the 98% confidence /2.0 a level) a model planet 
population using a planet distribution where N(a) oc constant out to a distance 
of 45 AU (further model assumptions discussed within). 

Subject headings: planets: extrasolar — instrumentation: adaptive optics — 
binaries: general 



1. Introduction 



While over 200 extrasolar planets have been detectecl^ over the last 11 years (mostly 
via the radial velocity technique), very few extrasolar planet candidates have been im- 
aged directly (for instance, 2 MASS 1207b (~8±3 Mj„p), Oph 1622B (~13±5 Mj„,J and 
CHXR 73 B (~12.5±8 Mj„p) Ichauvin et al]l2005al : Iciose et al.l[2007al : lluhman et alJbood : 
Brandeker etaDl2006h . The few candidates discovered of "planetary mass" < 13 Mj„p are 



companions to brown dwarfs and possess properties more similar to young brown dwarfs 
(separations > 50 AU; surface gravity g > 0.3) than to giant extrasolar planets orbiting 
sun-like stars. Based on their large (>50 AU) separations, these objects appear to have 
formed via a fragmentation process, more similar to brown dwarfs. Hence, to date no true 
images of extrasolar planets have been obtained. 



^http://exoplanet.eu/catalog.php, maintained by Jean Schneider 
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Theoretically, a large telescope (D > 6 meters) plus an adaptive optics (AO) system 
should be able to reach the photon-noise limit at 1" separations from the star with an 
hour of exposure time and thus attain the very high (>10^) contrasts necessary to image 
a young extrasolar giant planet. Thus, numerous adap tive optics surveys to directly detect 



extra solar planets have been completed (for instance, iKaisler et al.l l2003l : iMasciadri et al. 



20051 ). These surveys have yielded interesting contrast limits but no true extrasolar giant 
planet candidates. 

The difficulty in directly imaging extrasolar giant planets can be attributed to the unfor- 
tunate fact that bright quasi-static speckles (also known as super speckles) caused by slowly 
evolving instrumental aberr ations remain in ada ptive optics images even after adaptive optics 
correction (see for example iRacine et al.lll999l ). These super speckles evolve stochastically 
on relatively long (minute) timescales and also vary somewhat chro matically, producin g cor- 
related speckle noise which is very difficult to calibrate and remove (IRacine et al.lll999l ). For 
photon-noise limited data, the signal to noise S/N increases as t°'^, where t is the exposure 
time. Approximately speaking, for speckle-noise limited data, the S/N does not increase 
with time past a specific speckle-noise floor (limiting AO contrasts often to ~10^ at 0.5", 
Racine et al. 1999; Masciadri et al. 2005). More exactly, S/N does continue to increase with 
time, but as the speckle noise in successive frames becomes correlated, the N gain becomes 
considerably slower. Effectively independent exposures then have durations of many min- 
utes rather than a small fraction of a second (Racine et al. 1999). This correlated speckle 
noise is considerably above the photon noise limit and makes planet detection very difficult. 
Interestingly, space telescopes such as HST also suffer from limiting correlated speckle noise 
due to temperature variations which induce changes in the PSF (known as "breathing", 
Schneider et al.ll2003l ). 



Many observatories, including Gemini, Subaru, and the VLT, are currently building 
dedic ated planet-finding AO / coronagraph cameras in order to overcom e this speckle noise 
floor (jPohlen et al.l l2006l : [Macintosh et al.l l2006l : iTamura &: Lyd |2006| ) . A number of in- 
strument al speckle-attenuat i on methods have been propos ed, such as spectral differential 



(Marois et a. 



Thatte et al 



i maging (IRacine et al.l ll999l: iMarois et al.ll2000l . [2002 



20051). azimuthal differential iraaging 



20061) , integral fleld spectroscopy (ISparks fc Ford! l2002l : iBerton et al.l 12006 
20071 ). precise wavelength control methods such as those develop ed at the High 



Contrast Imaging Testbed (jTrauger et al.ll2004j ). focal plan e wavefront sensin g (jCodona fc Angel 
20041 : iKenworthv et al.ll2006h . and nulling interferometry (ILiu et al. II2006I ). 



The Simultaneous Differential Imagers at t 



our team ( Lenzen et~alll2004l . [2005l :[c 



attenuation technique (pioneered by 



ose et al 



le VLT and MMT, built and commisioned by 



Racine et al 



2005af). utilizes a spectral differential spec kle- 
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exploits a methane absorption feature at 1.62 /im (see Fig. [J) which is robustly observed 



in su bstellar objects with spectral type later than T3.5 (iGeballe et al.ll2002l : iBurrows et al. 
200ll ). SDI utilizes specialized hardware to image simultaneously inside and outside this 



methane feature with custom 25 nm filters (see Fig.[T]). Since the super-speckles are coherent 
with the starlight and both starlight and speckles have a fiat spectrum (see Fig. [1]) in this 
narrow wavelength band {5X / A ~1.6%), subtracting the "on" and "off" methane absorption 
images removes the starlight and its speckles, while preserving light from any substellar 
methane companion to the star. 

We have completed a 54 star survey with the SDI device at the VLT and MMT. Survey 
stars were chosen primarily according to proximity to the Sun (<50 pc) and youth (<300 
Myr, typically <100 Myr). We observed 47 young (<250 Myr) stars, 3 nearby stars with 
known RV planets, and 4 very close (<20 pc) older solar analogues. We obtained contrasts of 
AH>10 mag (5 a) at 0.5" for 45% of target objects at the VLT and contrasts of AH>9 mag 
(5 a) at 0.5" for 80% of our targets. The VLT SDI device is fully commissioned and available 
to the community and the MMT SDI device is a PI instrument with the A RIES camera. 



In cq i itrast, the dedicated pl anet-finding instruments such as Sphere and GPI (iDohlen et al. 



20061 : [Macintosh et al.l |2006| ) being built at the VLT and Gemini will not see first light for 
several years. Thus, as a precursor to planet surveys with these dedicated planet finding 
cameras, the results from the SDI devices are especiall y timely and rel evant, particularly to 
inform the large Gemini NICI survey starting in 2007 (ILiu et al.ll2005l ). 



2. The Simultaneous Differential Imagers at the VLT and MMT 

The VLT Simultaneous Differential Imager (henceforth SDI) was built at the University 
of Arizona by L. Close and installed in a special f/40 camera relay for the VLT AO camera 
CONICA built by R. Lenzen at the Max Planck Institute for Astronomy, Heidelberg. These 
were both installed at the VLT in August 2003. The MMT SDI was also built at the 
University of Arizona. In February 2004, it was installed in the ARIES f/30 camera built by 
D. McCarthy. Both devices are available to the observing communities of their respective 
telescopes. 



2.1. Hardware Considerations 



The SDI device consists of a custom double WoUaston, which splits the incoming AO 
beam into 4 identical beams (utilizing calcite birefringence to minimize non-common path 
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error - adding only <10 nm rms of differential non-common path errors per the first few 
Zernikes modes - Lenzen et al. 2004a). Each beam then passes through a narrowband filter 
with a central wavelength either on or off methane absorption. Three different filters were 
used; all filters were placed in different quadrants on the same substrate. SDI filters for the 
VLT and MMT were manufactured by Barr Associates. Filter wavelengths were chosen on 
and off the methane absorption feature at 1.62 /im and were spaced closely (every 0.025 /im) 
in order to limit residuals due to speckle and calcite chromatism. We used four filters FX, 
F2, F3a, and F3b with measured cold central wavelengths Fl = 1.575 /im, F2=1.600 ^m, 
and F3a=F3b=1.625 fim. The filters are approximately 0.025 fim in bandwidth (1.6%). 
The SDI filter transmission curves overlaid on a theoretical young planet spectrum (private 
communication, D. Sudarsky) are presented in Fig. [H 



2.2. Discoveries with the SDI Cameras 



The SDI device has already produced a numbe r of important scien tific results: the dis- 
covery of the important calibrator object AB Dor C (IClose et al.ll2005bl ) which is the tightest 
(0.16") low mass (0.090±0.05 M0, ~100x fainter) companion detected by direct imaging, 
the m ost detailed methane surface maps of Titan from the pre-Cassin i era (jHartung et al. 
2004), the discovery of e Ind Ba and Bb, the nearest binary brown dwarf (iMcCaughrean et al. 
2004r). the discovery of SCR 1845-6357B, a very close (3.85 pc) T6 brown dwarf ( IBiller et al. 
2006bl ). and evidence of o rbital motion for Gl 86B, the first known white dwarf companion 
to an exoplanet host star (IMugrauer fc Neuhauserll2005l ). In fact, the SDI device discovered 
all known brown dwarfs within 5 pc of the Sun. It has also set the best upper limit on the 
luminosity of the older (~1 Gyr) extrasolar planet around eEri. 



2.3. Observational Techniques and Data Reduction 

To ensure the highest possible signal to noise ratio and to maximize SDI speckle atten- 
uation, a complex data acquisition procedure was followed for each star. For each object 
observed, we saturated the inner ~0.1" of the star, thus providing a wide dynamic range and 
contrast down into the halo. Base exposure times (DIT) range from 0.3 to 20 s (typically 
this was > 2s to allow Fowler sampling at the VLT), depending on the H magnitude of the 
observed star. A number of exposures (NDIT) with the base exposure time are then coadded 
in hardware to produce a standard ~2 minute long base datum. An example raw datum is 
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presented in Fig. [2]q 

Base datum are then taken at a grid of dither positions (4x0.5" spacings with the MMT, 
5x0.5" spacings with the VLT). This dither pattern is then repeated at typically 2 telescope 
"roll angles" (where a "roll angle" refers to a different field derotator position / position 
angle (henceforth PA) settings). A subtraction of data taken at different roll angles further 
attenuates super-speckle residuals (since the weak residual speckles after SDI subtraction are 
instrumental features in the SDI optics which do not shift with a change in roll angle) while 
producing a very important signature "jump" in position for any physical companion (since 
a physical companion will appear to shift by the roll angle difference between datasets). For 
a space telescope such as Hubble (where the entire telescope can be rolled), a companion 
detected at the 5cr level in two different roll angles would be detected at the 7a level (a S/N 
gam of ~ V2) across the entire dataset (assuming roughly Gaussian statistics). This method 
is somewhat less effective with ground based telescopes where field rotation is provided by 
the field derotator rather than rolling the entire telescope (thus, super speckles from the 
telescope optics can appear to rotate by the roll angle as well). Nonetheless, observing at 
two roll angles provides us with two independent detections of a substellar companion at 
different locations on the detector, thus allowing us to rule out a "false positive" detection 
at an extremely high level of confidence - indeed, the only 3 faint companions (eind Bb, 
SCR 1845-6357B, and AB Dor C) ever detected with >5cr using SDI in more than one roll 
angle have all proven to be real. A typical observing block at the VLT then consists of the 
following series of : 1) ~10 minute long dither pattern taken with a roll angle of degrees. 
2) ~10 minute long dither pattern taken with a roll angle of 33 degrees. 3) ~10 minute long 
dither pattern taken with a roll angle of 33 degrees. 4) ~10 minute long dither pattern taken 
with a roll angle of degrees. A custom template was developed at the VLT to automate 
this process in each OB. 

Each base datum was reduced using a custom IDL pipeline (described in detail in Biller 
et al. (2006a) and Biller et al. (2006c)). This pipehne performs sky-subtraction, flat-fielding, 
and bad pixel removal, extracts a square aperture around each separate filter image, scales 
the platescale of each filter image so that the speckles in each filter fall at the same radii 
despite chromatic differences, scales the flux in each image to remove any quantum efficiency 
differences between the images, and filters out very low (>15 pixels) spatial frequencies by 
unsharp masking each image. Each filter image is then initially aligned to a reference image 



•^As with all our survey data, this was taken with the original SDI double WoUaston prism. In February 
2007, the original prism was replaced with a next generation prism which is cut in such a way that each 
subimage now subtends a whole quadrant of the detector chip. The new prism is also fabricated from YV04, 
a material which produces smaller chromatic errors at 1.6//m than the original calcite. 
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to within 0.25 pixels using a custom shift and subtract algorithm (Biller et al. (2006a,c)). 
One master reference image is used for each ~40 minute long dataset. After each of the 
filter images has been aligned to the reference image, we calculate 2 differences which are 
sensitive to substellar companions of spectral types T {Teff < 1200 K) and "Y" {Teff < 600 
K). The first is optimal for T spectral types: 

Difference! = Fl(1.575 /xm) - F3a(1.625 ^im) (1) 
The second is optimal for Y spectral types: 

Difference2 = F2(1.6 /im) - F3a(1.625 /xm) (2) 

An additional alignment is performed before the SDI subtraction; using the Fl image as 
our reference image, we align images Fl and F3a to within 0.05 pixels. A similar alignment 
is performed with images F2 and F3a, using the F2 image as the reference image. 

These differences are also somewhat sensitive to hotter substellar companions (L and 
early T spectral types), due to the fact that the platescale in each filter image has been 
scaled to a reference platescale to align the Airy patterns in each image. A real object (as 
opposed to a speckle) will not scale with the Airy pattern and thus, after scaling, will appear 
at a slightly different radius in each filter image. Subtracting images in different filters will 
then produce a characteristic dark-light radial pattern for a real object. This effect obviously 
scales with radius - at the VLT, an object at 0.5" will be offset by less than 1 pixel between 
filters, while an object at 1.5" will be offset by ~3 pixels, producing a very noticeable pattern. 
Thus, the SDI subtractions have a limited sensitivity to bright L and early T companions. 
We note that AB Dor C (AH ~ 5 mag) was detected at 0.15" (February 2004, Close et 
al. 2005) and 0.2" (September 2004, Nielsen et al. 2005) separations from AB Dor A even 
though AB Dor C has no methane absorption features (as is expected from its M5.5 spectral 
type. Close et al. 2007b.) 

We additionally calculate one further non-differenced combination sensitive to M, L, 
and early T companions: 

Broadband = Fl(1.575/im) + F2(1.6/im) + F3(1.625/xm) (3) 

After each datum is pipehned the data are further processed in IRAF. For each ~10 
minute long dither pattern, all three combinations described above and the four reduced 
filter images are median combined. Each 10 minute dataset is then differenced with the 
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following 10 minute dataset (taken at a different position angle). All roll-angle differenced 
images for each target object observation are then median combined to produce the final 
data product. 

A fully reduced ~30 minute dataset of AB Dor A (70 Myr KIV star at a distance of 
14.98 pc, V=6.88) from the VLT SDI device is presented in Fig. |3l Simulated planets have 
been added at separations of 0.55, 0.85, and 1.35" from the primary, with AFl(1.575/im) = 
10 mag (attenuation in magnitudes in the 1.575 /im Fl filter) fainter than the primary. For 
details and further discussion of these planet simulations see Section 3.4. 



3. The SDI Survey 

3.1. Survey Design / Target Selection 

Survey objects were selected primarily on the basis of youth and proximity. With a 
number of exceptions, our 54 survey objects are within 50 pc of the Sun and less than 250 
Myr in age. (The 9 exceptions include three somewhat older stars with known radial velocity 
planets, 2 more distant (<150 pc) stars with extreme youth indicators, and 4 older nearby 
young solar analogues which were initially misclassified as young objects.) Distances were 
obtained for 48 of our objects from Hipparcos parallax na easurements (parallaxes of >0.02", 



corresponding to distances <50 pc. lPerryman et al.lll997l ). Stars were age-selected according 



to two methods: 1) if possible, according to young cluster membership (and adopting the 
established age for that cluster) for clusters with well established ages such as the Beta 
Pic, TW Hya, AB Dor and Tuc-Hor moving groups or 2) according to other age indicators 
including the strength of spectral age indicators (for instance, the Li 6707, the Calcium H and 
K lines, and Ha emission) as well as from X-ray emission, variability, and rotational speed. 
As moving group ages are generally more robust than measurements for individual stars, 
we expect the ages of stars in these associations, on average, to have greater accuracy. Our 
survey covers stars in the Beta Pic, TW Hya, AB Dor, IC 2391, and Tucanae/Horologium 
moving groups. 

We select targets stars based on two overlapping criteria: 1) stars within 25 pc and 
younger than 250 Myr, and 2) stars within 50 pc and younger than 40 Myr (see Fig. H]). Our 
original list has been modified according to the amount of allocated time at the telescope, the 
unavailability of GTO targets, as well as severe weather constraints for the MMT portion of 
our survey. At the VLT, our observing runs spanned the months of August through February 
over 2004 and 2005. Thus, due to the spacing of observing runs, in the south, the survey is 
close to complete from ~17 - ~13 hours RA. At the MMT, we had two observing runs, one 
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Fig. 1. — SDI filter transmission curves overlaid on the theoretical spectrum (private com- 
munication, D. Sudarsky) of a young extrasolar planet (30 Myr, 3 Mj^p). Filters 1 and 2 
sample off the 1.62 /im CH4 absorption feature, while filter 3 samples within the absorption 
feature. In contrast, the spectrum of the K2V star eEri (Meyer et al. 1998) is flat across the 
whole wavelength band. Subtracting images taken in fllters "on" and "off" the methane ab- 
sorption feature will remove the star and speckle noise (which is coherent with the starlight) 
while preserving any light from giant planet companions. (Details of the complex SDI data 
pipeline are provided in Section 2.3.) 
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Fig. 2— Two minutes of raw SDI data from NACO SDI's 1024x1024 Aladdin array in 
the VLT CONICA AO camera (Lenzen et al. 2004). A number of electronic ghosts are 
apparent outside the four square filter apertures (each aperture is rotated by 30°); indeed, 
filter apertures were specifically selected to exclude these ghosts. Note that this is an image 
of the original Alladin array; the current SDI array has far fewer bad pixels. 
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Fig. 3. — Left: A complete reduced dataset (28 minutes of data at a series of rotator 
angles ("roll angles") - 0°, 33°, 33°, 0°) from the VLT SDI device. Simulated planets have 
been added at separations of 0.55, 0.85, and 1.35" from the primary, with AFl(1.575/im) 
= 10 mag (star-planet contrast in magnitudes) fainter than the primary. These planets 
are scaled from unsaturated images of the example star (AB Dor A) taken right before the 
example dataset (and have fluxes and photon noise in each fllter appropriate for a T6 effective 
temperature). Past 0.7", the simulated planets are detected in both roll angles with S/N > 
10. Observing at two different roll angles produces two independent detections, and hence 
makes the chance of detecting a "false positive" almost null. Right: Standard AO data 
reduction of the same dataset. Filter images have been coadded (rather than subtracted), 
flat-flelded, sky-subtracted, and unsharp- masked. Simulated planets have been added with 
the same properties and at the same separations as before. None of the simulated planets 
are clearly detected in the standard AO reduction. Additionally, many more bright super 
speckles remain in the fleld. 
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in May 2005 and one in February 2006. Thus, in the north, the survey is complete for the 
RA range 11-21 hours. 

Survey objects are presented in Table 1. A detailed table of observations is presented 
in Table 2. Survey objects are plotted as a function of distance and age in Fig. HI Our 
"median" survey object is a K star with an age of 30 Myr and at a distance of 25 pc. 



3.2. The Performance of the SDI Filters as Spectral Indices 



It is important to carefully consider the expected strength of the 1.62 fim methane 
absorption break utilized by the SDI device. The stronger the break strength, the more 
companion light is preserved after SDI filter subtraction. For a candidate object with a 
weak break strength, SDI subtraction may effectively attenuate the candidate object itself, 
rendering it undetectable (although, at separations > 0.15", a bright object may still be 
detectable due to the characteristic dark-light radial pattern produced by any real object 
after pipelining, see Section 2.2.) 

To determine the methane break strength expected for a candidate object (and thus, 
the expected performance of SDI for that candidate), we define an SDI methane spec- 
tral index calculated from our SDI Fl(1.575 /xm) and F3( 1.625 /zm) filter images (simi- 
lar to the methane spectral index defined by iGeballe et al.l 120021 ). 



index{ 



Fl. 
F3 



A2=1.5875^im 
-1.5625fim 



SxFl{X)d\ 



A4=1.6125/im 
Xi=1.6375fj.m 



SxF3{X)dX 



(4) 



Each SDI filter was manufactured by Barr Associates to have a precise bandwidth of 
0.025 yum, so the wavelength intervals (A2 - Ai = AA = A4 - A3) in the numerator and 
denominator have the same length for the SDI methane index. 

We calculated SDI spe ctral indices for the fo ur brown dwarfs which have been observed 
with SDI - the T6 Gl 229B jNakaiima et al.lll995h . the T5.5 SCR 1845B jBiller et al.ll2006bh 



and eind Ba-Bb (T6 + Tl) jMcCauehrean et al 



20041 ). Since we only possess SDI data on 



a limited number of T dwarf s, we calculated th e same SDI spectral indices from spectra of 
56 L dwarfs and 35 T dwarfs (iKnapp et al. I I2OO4I ) in order to evaluate the performance of the 



SDI for a wide range of L and T dwarf objects. Spectra for these objects were obtained from 
Sandy Leggett's L and T dwarf archive Q In order to make an accurate comparison, SDI 



^http://www .jach.hawaii.edu/~skl/LTdata. html 
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Fig. 4. — Age vs. distance for our survey stars. Spectral types are delineated by plot symbols. 
Objects were selected according to youth and proximity to the Sun. 45 of our survey objects 
are within 50 pc of the Sun and less than 250 Myr in age. Of the remaining objects, 2 are 
very young (<10 Myr), somewhat more distant (<150 pc) objects, 3 are nearby stars with 
known RV planets, and 4 are nearby solar analogues (<20 pc) that were initially misclassified 
as young. We selected targets according to two overlapping criteria (shown on plot as solid 
black lines) 1) stars within 25 pc and younger than 250 Myr and 2) stars within 50 pc and 
younger than 40 Myr. Stars were age-selected according to association membership, or, in 
the case of unassociated stars, age indicators such as the strength of the Li 6707 A line. 
Calcium H and K lines, H a emission. X-ray emission, etc. Distances were obtained from 
Hipparcos parallax measurements (parallaxes of >0.02"). Our "median" survey object is a 
K star with an age of 30 Myr and at a distance of 25 pc. 



-14- 



filter transmission curves were convolved into these calculations (see Fig. [T]). Since we have 
full spectr al data for these obj ects, we also calculated the 1.62 fim methane spectral index 
defined by lGeballe et al.l (120021 ) . which were found to be similar to our SDI methane spectral 
indices. SDI methane spectral indices are plotted for both the M9 and T6 components of 
SCR 18 45, the T dw a rfs G l 229B, eind Ba, eind Bb, and 94 other L and T dwarfs in 
Fig. [5l iGeballe et al.l (120021 ) note that Gl 229B has an anomalously high methane index 
for its spectral type and assign a large uncertainty to Gl 229B's spectral type - T6±l - 
which is also refiected in its anomalously large SDI spectral index compared to other T6 
dwarfs. From this analysis, we conclude that the SDI device can effectively detect objects 
with spectral type later than T3. Since T dwarfs with spectral type earlier than T3 are 
relatively uncommon compared to later T dwarfs, the SDI device can effectively detect the 
full range of extrasolar giant planet / brown dwarf spectral types of interest. According to 
the models of Burrows et al. 2003 and Marley et al. 2006, planets >10 Myr old should 
possess Te// < 800 K and have spectral type of T8 or greater. 



3.3. Contrast Limits and Minimum Detectable Planet Separation 

To determine the range of possible star-planet contrasts achieved in our survey, we 
generated noise curves as a function of radius for every survey star. We tested three different 
methods of generating noise curves: 1) translating a 6x6 pixel (0.1"x0.1") box along a 
particular radial trajectory away from the center of the star image (typical PSF FWHM 
was 3-5 pixels) then calculating the standard deviation in the box at each point along this 
trajectory, 2) averaging noise curves generated along four such trajectories, and 3) calculating 
the standard deviation within annular regions 6 pixels in width centered on the primary 
PSF (spider diffraction spikes were not masked out in this case because they are already 
well removed by the spectral difference). Noise curves generated in these three manners are 
presented for a set of 6 typical program stars (AB Dor, DX Leo, GJ 182, AB Pic, GJ 799A, 
and GJ 799B) in Fig. [HI In general, all three methods produce remarkably similar noise curves 
and are equally suitable for characterizing the noise properties of an observation. However, 
we choose to utilize the single trajectory method because it best simulates the particular 
signal to noise issues encountered when searching for faint companions among super-speckles 
of similar intensity and FWHM (since it preserves pixel to pixel noise variations due to 
super-speckles). The annular method averages out speckle noise properties azimuthally. 
This produces somewhat unrealistic results in the case of a faint companion search where 
one is concerned only with the speckle structure within the local area of a candidate faint 
companion - speckle structure on the other side of the image is unimportant. In addition, 
we have tried to choose very "typical" trajectory per star - ideally, trajectory to trajectory 
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Fig. 5. — SDI methane spectral indices for the T dwarfs SCR 1845B, Gl 229B, eind Ba, and 
eind Bb (from Biller et al. 2006b). As a comparison, S DI methane spectral indices calculated 
from spectra for 94 L and T dwarfs (spectra from iKnapp et al.l |200J) are overplotted. 
SCR 1845B, Gl 229B, and eind Bb show strong methane indices, whereas eind Bb (Tl) 
is relatively c o nstari t in flux across the SDI filters and has a much lower methane index. 
Geballe et al.l (j2002[ ) note that Gl 229B has an anomalously high methane index for its 
spectral type. While iGeballe et al.l (120021 ) find an overall spectral type of T6±l for Gl 229B, 
they assign Gl 229B a spectral type of T7 based on the methane index (which we adopt 
here) . 
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variations will average out across the entire survey. 

Noise curves for each program star were calculated along a trajectory 45° from the image 
X axis in the first quadrant. The 45° was selected as one of many possible representative 
trajectories which was unaffected by instrumental effects such as spider arms, vibrations 
along azimuth or altitude mounts, etc. At each point along this trajectory, the standard 
deviation was calculated (except for the PSF noise curve, for which the mean was calculated). 

A fully labeled example noise curve for the star DX Leo is presented in Fig. [71 Noise 
curves were generated for a number of cases for each object. First, a noise curve was generated 
for the full reduced and differenced SDI data (labeled SDI data curve) (Fl (1.575 /im) - 
F3a(1.625 /im) for two roll angles). A PSF noise curve curve was generated from a median 
combination of all the Fl (1.575 /im) filter images for each dataset weighted according to 
the number of exposures, dithers, and roll angles in the dataset. To recreate the equivalent 
observation without using the SDI technique (and thus characterize the performance of 
SDI compared to conventional AO techniques), an "optimized conventional AO" curve was 
generated by combining images from all three filters at each roll angle: 



Broadband = Fl(1.575/xm) + F2(1.6/im) + F3(1.625/im) (5) 

then unsharp masking to remove low spatial frequencies, and subtracting the "Broad- 
band" combinations at different roll angles from each other. 

To characterize the noise level in each observation, we calculated an SDI noise curve, 
which is a combination of photon-noise, flat-field noise, and read noise. Per exposure: 

'^SDI = \j (^photon + '^}lat field + ^readnoise (6) 

Photon-noise was calculated as: 



'^photon yf^ electrons \' ) 

Readout noise for the CONICA detector at the VLT in Fowler sampling mode is 1.3 ADU 
(analog-to-digital unit). The gain for the latest CONICA detector in the Fowler sampling 
mode is 12.1 electrons/ ADU so areadnoise = 15.73 electrons. 

NACO and ARIES flat fields were found to be accurate to about 1%, so fiat-field noise 
was estimated as: 



^ flatfield ^^electrons 



(8) 



where e=0.01. The total noise for a full observation (4-5 dithers, 2-4 roll angles) was 
then calculated by weighting the SDI noise per exposure by the number of exposures (NDIT 
X number of dithers x number of roll angles): 

o'sDi_fuiiobs = crsDi\/ NDIT X (number of dithers) x [number of roll angles) (9) 
The PSF curve for a full observation was similarly weighted: 

PSF = (medianPSF) x NDIT x [number of dithers) x [number of roll angles) (10) 

For the sample curve shown in Fig. [71 the SDI data is "flat-field" limited within 0.5" 
of the star. From 0.5" onwards, the SDI data is photon noise limited, approaching the 
read-noise limit at separations > 2" . 

We converted our noise in electrons to attainable contrasts in magnitudes in the Fl (1.625 
/im) filter - contrast plots in Amag are presented for all non-binary survey objects in Figs. [S] 
to [m according to the H magnitude of the primary for the VLT and according to observing 
run for the MMT. For every observation which possesses an unsaturated acquisition image 
(typically 10x0.1 s images taken over ~30 s), the stellar peak in the unsaturated acquisition 
image was used to scale the saturated stellar peak in the saturated data images and thus 
attain accurate contrasts in magnitudes. For observations lacking an unsaturated acquisition 
image, contrast curves for other stars which had similar peaks, read noise values, and shape 
to the contrast curve in question were selected from the library of contrast plots in electron 
units. The peaks utilized for these matching contrast curves were then used to scale the 
observation missing an acquisition image. A peak of 2.2x10^ was adopted for eEri (Kellner 
et al. 2007, Janson et al. 2007) and eind A (GeiBler et al. 2007). We present contrast curves 
for 48 stars in this paper; the remaining six survey stars were either very close binaries, 
making it difficult to generate a contrast curve, or had particularly low quality datasets. 

For the VLT data, attainable contrast depends on primary star H magnitude as well 
as seeing FWHM and Strehl ratio during the observation. For the brightest stars in the 
survey (II<4.5), we attain 5o" contrasts of AF1~12 mag at separations of >1" from the star. 
For the faintest survey stars, we only attain 5(t contrasts of AF1~10 mag >1" from the 
star. However, considerable spread in attained contrast is observed in each H magnitude 
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Fig. 6. — Comparison of Noise Curves generated in 3 different manners for a set of 6 typical 
program stars (upper left: AB Dor, upper right: DX Leo, middle left: GJ 182, middle 
right: AB Pic, lower left: GJ 799A, lower right: GJ 799B). Noise curves were generated 
by: 1) translating a 6x6 pixel (0.1"x0.1") box along a particular radial trajectory away 
from the center of the star image (typical PSF FWHM was 3-5 pixels) then calculating the 
standard deviation in the box at each point along this trajectory, 2) averaging noise curves 
generated along four such trajectories, and 3) calculating the standard deviation within 
annular regions 6 pixels in width centered on the primary PSF (spider diffraction spikes 
were not masked out in this case because they are already well removed by the spectral 
difference). In general, all three methods produce remarkably similar noise curves and are 
equally suitable for characterizing the noise properties of an observation. Since it preserves 
pixel to pixel contrast variations due to speckle noise, the single trajectory method better 
simulates the S/N issues encountered in searching for faint companions. 
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Fig. 7.— Sensitivity curve for DX Leo (18 pc, KOV, 115 Myr, V=7.05, H=5.242). This 
is 28 minutes of VLT SDl data. The CONICA PSF curve is the median combination of 
all the Fl (1.575 ^m) filter images for this dataset (with a gain correction apphed which 
accounted for the number of exposures, dithers, and roll angles). The "optimized conven- 
tional AO" curve was generated by averaging images from all three filters at each roll angle, 
unsharp masking to remove low spatial frequencies, then subtracting the combinations at 
different roll angles from each other. The "measured SDl" data curve is the full reduced 
and differenced SDl data for this object (Fl(1.575 /im) - F3a(1.625 /im) for two roll an- 
gles). The "theoretical SDl noise" curve is calculated from photon noise (long dashed green 
curve), fiat-field noise (short dashed black curve), and read noise (solid black line) added 
in quadrature. Within 0.5", the SDl data is "fiat-field" noise limited. (In reality, we are 
limited by super speckle residuals within this radius. Our fiat fields are accurate to the 
~1% level, but the speckle residuals <0.5" vary more than this and thus dominate the SDl 
noise.) From 0.5" onwards, the SDl data is photon-noise limited, asymptotically approach- 
ing the read-noise limit at separations > 2". For a complete set of sensitivity curves, see: 
http: //exoplanet. as. arizona.edu/~lclose/SDI.html 
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bin - most likely due to variations in observing conditions (seeing, Strehl ratio, etc.) across 
multiple observations. To quantify the effect of seeing on attainable contrast, in Fig. [15] we 
plot the seeing FWHM (averaged over the observation - the error bars on seeing are the 
seeing variations as measured by the standard deviation of the seeing over each observation) 
vs. attained 5a contrast at 0.5" for 10 of the stars presented in Fig. [9] with H magnitudes 
between 4.5 - 5.5. For this sample of stars with similar H magnitudes, achievable contrast is 
roughly inversely proportional to the seeing FWHM. A fair amount of scatter is apparent in 
this plot and is due in part to seeing variations over the course of each observations. Seeing 
FWHM can vary considerably over the 20-40 minute timescale of a typical SDI observation, 
affecting the AO system performance and thus the achievable contrast. 

However, higher attained contrast does not necessarily translate across the board to a 
lower minimum detectable planet mass. Although one might be able to attain a very high 
contrast {5a contrast >11 mag at 1" limited by photon noise) for a bright young A star, 
one would have more luck searching for low luminosity planets around an intrinsically faint 
young M star {5a contrast ~9 mag at 1" limited by read noise), since the inherent contrast 
difference expected between star and planet is considerably smaller. We obtained contrasts 
of AH>10 mag (5 a) at 0.5" for 45% of target objects at the VLT and contrasts of AH>9 
mag (5 a) at 0.5" for 80% of our targets. This is more a statement on the spectral types in 
our sample than a performance related issue. 

In general, the MMT SDI device performed at a slightly lower level than the VLT SDI 
device - attaining 5a contrasts 0.5-1 magnitude less than those achieved at the VLT for 
similar separation and primary star H magnitude. The lesser performance of the MMT 
system can be attributed to two factors. First, the diameter of the MMT is 6.5m versus 
the VLT which has an 8.2 m diameter - resulting in a considerable decrease in sensitivity. 
Additionally, the seeing sampled by the MMTAO system was not as stable as for the NACO 
AO system - Strehl ratios often changed dramatically over an observation, limiting the 
attainable contrast. However, the MMT SDI results still probe a higher contrast regime at 
separations <1" than is possible with standard AO techniques. 

In order to determine what objects realistically can be detected for our survey stars, 
we must convert between our instrumental Fl(1.625 /im) filter magnitudes and H band 
magnitudes and then compare the H magnitudes to those expected from models of young 
planets (such as Burrows et al. 2003). To accomplish this, the spectra of both the primary 
and secondary components of each target must be taken into account. To convert from 
our Fl filter magnitudes into calibrated H band magnitudes we must calculate the H band 
magnitude offsets for both the primary star and a potential methane companion (Offset^ 
and Offset^ respectively): 
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AH = Ha-Hb = {OffsetB+FlB)-{OffsetA+FU) = [O f f sets-O f f setA)+AFl (11) 

For primary stars with spectral types F-K, we assume that the star has very httle 
chromatic variation within the middle of the H band, so Offset^ is zero (see Fig. [1]). For 
lower mass M stars, which are very red, the magnitude offset is not negligible. To take an 
extreme example, a very low mass M8 primary will have a magnitude offset of Offset^=- 
0.12±0.08 mag (calculated using the spectrum of the M8 star VBIO, an H transmission 
curve, and our Fl filter transmission curve). The latest stars in our survey have spectral 
type MO- M5, so Offsetyi will be <0.1 mag for these cases. 

Any T3 or later companion to one of our survey stars will be blue compared to the 
primary and will appear "brighter" in the Fl filter than in the H band (in other words, it 
will have a higher "flux" in the Fl filter (7^ photons per unit bandwidth) - see Fig. [1]) - so 
Offsets will definitely be non- n egligible. We calcu lated Offset^ for 18 objects with spectral 



types of T4.5-T8 (spectra from lKnapp et al.ll2004l ). then averaged together by spectral type 



to derive an average offset for each spectral type. For a T5 companion, Offset-rs = 0.5±0.05 
mag, for a T6 companion, Offset^e = 0.6±0.07 mag, and for a T8 companion, Offset^g 
= 0.87±0.04 mag. While we do not convert our full AFl contrast plots to AH contrast 
plots, for every survey star we calculate limiting AH contrasts (5cr values), at 0.5" and 1.0", 
equivalent separation in AU, apparent H magnitude, and absolute H magnitude for a T8 
spectral type companion (since extrasolar planets are expected to have spectral type > T8, 
Burrows et al. 2003). These results are presented in Tables 4 and 5. However, it is difficult 
to translate our absolute H magnitudes into model planet masses since we have assumed a 
T8 spectral type in our conversion between AFl and AH contrasts - but a companion which 
actually has the limiting absolute H magnitude we find (combined with the known age and 
distance of the system) may have a very different spectral type. 

Since we cannot translate our H magnitudes directly into planetary mass companions, 
we followed the analysis of Masciadri et al. (2005) and translated theoretical planet models 
(Burrows et al. 2003, Baraffe et al. 2003) into H magnitudes then determined the minimum 
separation at which such a companion could be detected (at the 5cr level) in our survey. 
The minimum separation at which a 5 Mj„p or a 10 Mj^p companion could be detected for 
each of our survey stars is shown in Table 6. Using the Burrows et al. (2003) models, for 
50% of the cases in our survey we detect no 5 Mj„p planets at separations larger than 18.6 
AU and no 10 Mj„p planets are separations larger than 7.5 AU. While these numbers are 
comparable to those found in Masciadri et al. (2005), our current survey actually attains 
higher contrasts on a case by case basis than Masciadri et al. (2005). Our median survey 
object has an age of 50 Myr whereas the median survey object of Masciadri et al. (2005) 
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has a considerably younger age of 12 Myr - the star-planet contrast is less at younger ages, 
thus one would expect a younger object to have a lower minimum separation at a given 
attained contrast than a similar but older object. For the 10 objects in common between 
the surveys, our survey attains lower minimum separations for 8 out of 10 objects (we note 
also that the two objects for which we did not attain lower separations were particularly 
low quality SDI datasets). Minimum detectable separations for a 5 IsAjup object for the 10 
objects in common are plotted in Fig. [16] (using the ages adopted by Masciadri et al. 2005). 
Our survey is generally more sensitive than Masciadri et al. (2005) on shared stars because 
the SDI technique allows us to achieve higher contrasts closer to the star (separations of 0.3" 
- 1.0") compared to the deep broad-band imaging technique of Masciadri et al. (2005), thus 
allowing us to potentially detect companions at tighter separations. We also shared 4 survey 
objects in common with Lowrance et al. (2005) and 1 object (eEri) in common with Luhman 
and Jayawardhana (2002). In all of these cases, our limiting contrasts at 0.5" (AH~10-11 
mag) are considerably higher than those attained in these previous surveys (AH~6.5-7.6 
mag), thus we are sensitive to planets at much smaller separations with SDI. 



3.4. Survey Completeness 

One would not expect a planet to be detectable at all phases of its orbit - to really 
understand the types of planets to which we are sensitive, we must take orbital motion into 
account and translate separations on the sky into orbital semi-major axes (a). To this end, 
we generated contour plots of fractional completeness as a function of mass and semi-major 
axis. For every survey star, we simulate 10000 planets for each combination of mass and 
semi-major axis. Eccentricities are drawn from a distribution of eccentricities consistent 
with known radial velocity planets. Standard distributions were used to randomly compute 
viewing angle and orbital phase, giving an instantaneous separation between star and planet. 
We use the distance, age, spectral type, and H-band magnitude of the star, and luminosity 
as a function of mass, calculated from the Burrows et al. (2003) models, to provide each 
simulated planet a separation on the sky in arcseconds, and an H-band flux ratio compared 
to its parent star. Combining this with the SDI contrast curve for each star in the survey, we 
can then determine the percentage of simulated planets detected as a function of mass and 
semi-major axis for each survey star. Contour plots for a set of 4 typical program stars (AB 
Dor, DX Leo, GJ 182, and GJ 799B) are presented in Fig. [T71 Note that we conservatively 
assume only T-type objects can be detected, hence masses > 10 Mj„p are not considered for 
many young targets. The value attached to each contour level defines the completeness of 
our observation to detecting (at the 5cr confidence level) a planet with the given semi-major 
axis and mass. It is worth noting that the only assumptions necessary for the generation of 
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these plots is the eccentricity distribution of planets and the Burrows et al. 2003 models. 



We use this method to summarize our survey completeness in Fig.[TSl Having computed 
the completeness for each star to planets at various masses and semi-major axes, we take 
slices at representative values of the semi-major axis, and present the number of stars in 
our 54 star survey which are at least 50% complete to such a planet. Our survey places the 
strongest constraints on planets between 6-10 Mj^p with semi-major axes between 20-40 AU. 
With 20 such stars (with 50% or greater completeness in this mass/semi-major axis range) 
surveyed without a detection of a planet, a simple way of interpreting our results (though 
without statistical rigor) is that we would expect the frequency of such planets to be of order 
10% or less. 



The evolutionary models of iBurrows et al.l (120031 ) utilize a "hot start" initial condition 
which, while appropriate for b rown dwarfs, is possib ly significantly different from the actual 
initial origins of planets. The iBurrows et al.l (120031 ) models begin with a high-temperature, 
high-entropy hydrogen-helium sphere which is allowed to radiate and cool over time. In 
contrast, a planet forms when gas accretes onto a rocky core, according to the core-accretion 
models of Id a and Lin (2005) and the disk instability models of Boss (2003). Recently, 
Mar ley et al.l (120061 ) simulated model planets with more realistic (lower entropy) initial con- 
ditions. These model planets have significantly lower luminosities at young ages (<1 Gyr). 
Model planets also converge to the "hot start" evolutionary tracks at different times accord- 
ing to mass - a 1 Mj^p model converges to traditional tracks by 20 Myr, while a 10 Mj^p 
requires up to 1 Gyr to match traditional tracks. Currently, H band magnitudes for these 
models are not yet available, but will be available in Spring 2007 (private communication, J. 
Fortney). When H band magnitudes are available, we will repeat this analysis using these 
new models. 



3.5. Sensitivity Case Study: AB Dor with Simulated Planets 



Since our survey data are highly saturated in the core of the image, it is difficult to 
place simulated objects in our data with a high degree of positional accuracy, as there is 
no external reference for position between data taken at different dithers and roll angles. 
However, as part of the SDI survey observations, our team di scovered a close-in (0.156") 
companion (hereafter AB Dor C) to the young star AB Dor Jciose et all l2005bh. While 



this co mpanion is a very low mass M star (0.090±0.005 M^^^, M5.5±l. IClose et al.ll2005bl . 
2007bl ) and hence, does not possess methane absorption features, it it still clearly detected 
in our SDI data. In our second AB Dor dataset where AB Dor C is separated from its 
primary by 0.2" (Nielsen et al. 2005), the AB Dor C source can be used to our advantage 
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as a reference position from which to offset - allowing us to add simulated planets into 
this dataset with highly accurate positions and relative fluxes independent of our "pipeline" 
calculated centroids. 

Simulated planets were produced by scaling ~10x0.1 s unsaturated images of AB Dor 
A taken right before the example dataset. Planets were simulated with AFl(1.575yum) = 
9, 10, 11, and 12 mag and with methane break strengths appropriate for T5, T6, and T8 
spectral types. Methane break strengths were calculated using the methane spectral index 
defined in Section 3.2. Photon noise and zero points appropriate for each object was added 
using the IRAF artdata/mkobject tool. The photometric zero point was calculated from AB 
Dor C. 

A fully reduced 28 minute dataset of AB Dor A (70 Myr KIV at a distance of 14.98 pc, 
V=6.88) from the VLT SDI device is presented in Fig. [12] with simulated planets added at 
separations of 0.4", 0.6", 0.8", 1.0", 1.2", 1.4", 1.6", 1.8", 2.0, and 2.2" from the primary 
(AFl(1.575^m) = 9, 10, 11, and 12 mag and spectral type T8). Past 0.7", the AFl(1.575/im) 
= 10 simulated planets are detected with S/N > 10. The 2.2" object falls off the edge of the 
aperture in several dithers and thus appears somewhat attenuated compared to the other 
simulated objects. Maximum achievable companion contrast at the 5a level as a function 
of distance from the star is plotted in Fig. [201 The residual noise curve for this star (see 
section 3.3) is also overplotted. Contrast curves (Sex) calculated with both techniques agree 
well with each other. Using the magnitude offsets developed in section 3.4, we convert our 
AFl(1.575/im) contrasts into AH for each spectral type. We adopt Offset^ = mag. Offset^ 
= 0.5 mag for a T5 object. Offsets = 0.6 mag for a T6 object, and Offset^ = 0.87 mag for 
a T8 object. AH vs. separation in arcsec is presented in Fig. | 



AFl contrasts were translated into planet masses using the 100 Myr models of lBurrows et al 



(l2003l ). According to the 100 Myr old model, objects with mass < 10 Mj^p will have Teff 
< 900 K - these objects are reliably of spectral types later than T7 (temperature scale from 
Burgasser et al. 2006). Thus, we adopt the T8 spectral type curve for this analysis. AB 
Dor has a like ly age of 50-70 Myr ( Nielsen et al. 2005, Close et al. 2007b) - we interpolate 



the models of lBurrows et al.l (120031 ) to derive masses at these ages as well. The minimum 
detectable planet mass as a function of distance from the star is plotted in Fig. [251 Adopt- 
ing an age of 70 Myr for AB Dor A, we can detect a 5 Mjup planet 12 AU from the star. 
However, as noted above, the Burrows et al. 2003 models utilize a hot start initial condition 
which may be inappropriate for a young planet. The Marley et al. (2006) models utilize 
more appropriate initial conditions and when H band magnitudes become available for these 
models, we will repeat this analysis. 
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3.6. Comparison with Other Direct Detection Methods 

We believe that our SDI images are the highest contrast astronomical images ever made 
from ground or space for methane rich companions <1" from their star. To substantiate 
this claim, we compare our SDI contrast curves with those produced using a variety of 
other competing methods (Azimuthal Differential Imaging (ADI), Marois et al. 2006, Lyot 
Coronagraph, Hinkley et al. 2007, HST NICMOS, Schneider et al. 2003, K-band Keck AO, 
Schneider et al. 2003, and NACO deep imaging in the Ks band, Masciadri et al. 2005). 
Comparison contrast curves are presented in Fig. [23l Apart from the Lyot and NICMOS 
curves, all curves are from >8m class telescopes. For ease of comparison, we convert our 
AF1=1.575 /im SDI contrast curve into the equivalent AH contrast appropriate for a T8 
spectral type companion. For methanated companions, SDI provides improved contrast by 
1-4 mag within 1" as compared to other methods. 

3.7. New and Confirmed Close Binary Stars 

A number of close binary stars were discovered or confirmed during our survey. In 
Table 7, we present separations and position angles measured from unsaturated SDI images 
of these stars acquired before each full SDI dataset was taken. These values are meant as 
estimates, hence, no error estimate is provided. We discovered close stellar companions to 
HIP 9141 (0.15" measured SDI separation), AB Dor A (0.16" measured SDI separation, 
see Close et al. 2005a), HD 48189A (0.14" measured SDI separation), HD 135363 (0.26" 
measured SDI separation) and CD-64 1208 (0.18" measured SDI separation). The <0.5" 
separation between the primary stars and these object makes it highly improbable that 
they are background objects. Additionally, we confirmed the close binary RXJ 1243.6-7834 
(0.068" measured SDI separation) discovered by Brandner et al. (2000), the visual double 
LH 98 062 (2.4" measured SDI separation) discovered by Mochnacki et al. (2002), the 
spectroscopic binary TWA 4 (0.78" measured SDI separation) discovered by Torres et al. 
(1995) and the close binary EK Dra (0.67" measured SDI separation) discovered by Metchev 
and Hillenbrand (2004). 

3.8. Candidate Identification / Elimination 

Survey data were examined for planet candidates by eye and also using automated 
detection algorithms; generally, the human eye proved more effective for detecting candidates. 
We identified 8 very tentative planet candidates at the VLT which passed the following tests: 
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1) Candidate must appear at the appropriate positions in the full reduced data. (i.e. 
candidate image position must jump by the appropriate roll angle.) 

2) Candidate must appear (at least marginally) at the appropriate position in each of 
the separate roll angle images 

3) Candidates detected in the Fl(1.575 /im) - F3a(1.625 /xm) difference should also be 
detected in the F2(1.6 /im) - F3a(1.625 /im) difference as well. 

These extremely tentative (<2cr) candidates are noted in the comments column of Table 
1, with the predicted mass (from the models of Burrows et al. 2003) and separation had it 
been real. No candidates were detected with > 3cr. None of the 8 tentative candidates were 
detected at a second epoch, thus the survey reached a null result for extrasolar planets at 
the ~3(T level and certainly at the 5 a level analyzed here. 



3.9. Planet Detectability 

To determine what sort of planets we can detect in this survey, we converted our contrast 
curves in Amag units into minimum detectable mass vs. separatio n (assuming a l a te T to 



early Y spectral type for all possible objects and using the models of iBurrows et al.l (120031 )). 
We calculated minimum detectable mass vs. separation for all stars with contrast curves 
in Figs. [8] to Figs. [TH minimum detectable mass vs. separation is presented for a set of 
four typical survey stars (AB Dor, DX Leo, GJ 182, and GJ 799B) in Fig. [211 However, 
to detect an object of any given mass requires that such an object exists around its parent 
object! The likelihood of detecting any object at a given radius is a combination of the 
minimum detectable mass for the parent star at that radius and the likelihood of such an 
object existing. Therefore it is very important to fully characterize and understand the 
expected distribution of objects around each survey star. The results of the survey then also 
constrain the possible distribution of extrasolar planets as a function of radius. 

To this end, we ran detailed Monte Carlo simulations to characterize the ensemble of 
planets expected to exist around each star. We conduct a similar simulation to that used 
to produce the contour plots of Fig. [T71 as described in Section 3.4 (these simulations are 
described in much more depth in Nielsen et al. 2006). In contrast to the production of the 
contour plots, we simulate 10^ planets instead of 10*^, and mass and semi-major axis are 
now assigned distributions of their own. The mass and semi-major axis distributions, like 
the distribution for eccentricity, are produced by considering the population of published 
radial velocity planets (e.g. Butler et al. 2006), with mass and eccentricity both chosen to 
fit the histograms from observed planets. Semi-major axis has been observed to follow a 
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distribution of N(a) oc for radial velocity planets (Wright et al. 2005). Since the radial 
velocity method has an inherent bias toward close-in planets (which have shorter orbital 
periods and larger radial velocity amplitudes), we attempt to correct for this by assuming 
a power-law distribution that is constant in semi-major axis - i.e. N(a) oc constant. We 
consider the results of Fischer and Valenti's (2005) volume-limited sample, and choose an 
outer limit for the semi-major axis distribution such that, for stars in the metallicity range 
in our sample, each star is expected to host one planet. This is done by integrating the 
semi-major axis distribution from 0.02 AU (corresponding to HD 41004Bb, the closest-in 
exoplanet known thus far) to 2.5 AU, the detection limit for the sample of Fischer and 
Valenti (2005), then noting the fraction of stars with planets in the metallicity range (-0.5 < 
[Fe/H] < 0.25) of our target stars (4.1%) and choosing an upper cut-off to the distribution 
when the integral reaches 100%. This gives us a constant probability distribution for semi- 
major axis between 0.02 and 45 AU that contains the same number of planets found in the 
<2.5 AU radial velocity survey. 

The ensemble of simulated planets is shown for our set of four typical stars in Fig. [2H 
Simulated planets which are detected are plotted as blue dots and those that remain unde- 
tected are plotted as red dots. In addition to the contrast plot, we also consider a planet 
"undetectable" when its apparent H magnitude drops below 21 mag (a limit set by our total 
integration time), or when the planet's temperature rises above 1400 K (given as a function 
of age and planet mass by Burrows et al. 2003). Above this temperature, the strength of the 
1.62 /im methane break weakens to the point that the SDI method loses effectiveness. Since 
we assume that each program star possesses exactly one planet that follows the distributions 
given above, we can assign a detection probability for that star from the percentage of the 
simulated planets that are detectable at the 5(T level. For our 48 program stars (consisting 
of 40 stars with ages >250 Myr and closer than 50 pc, 1 10 Myr old star at a distance of 
67 pc, 3 stars with known RV planets and 4 nearby solar analogues) which possess contrast 
curves, the average detection probability is 8.0%, the median detection probabihty is 4.1%, 
and the maximum detection probability is 47%. We have chosen to leave the older stars in 
this sample in our statistics even though their detection probabilities are essentially zero. 
Integrating over the probability distribution of our program stars, in Fig. [25] we plot the 
number of planets we expect to detect as a function of total stars observed, ordering the 
results so that the best stars (highest detection probabilities) are considered first. For the 48 
stars in our surveys for which we acquired contrast curves, we expect to detect a total of 3-4 
planets (3.8 to be exact) based on the above assumptions. Thus, our survey null detection 
rules out this exoplanet distribution at the 98% (2. Oct) level. It is important to note that 
this null result shows that this particular combination of assumptions (mass distribution, 
eccentricity distribution, constant semi-major axis distribution, upper limit to semi-major 
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axis at 45 AU, assumption that each star has a planet, and the mass-luminosity conversion 
from the models of Burrows et al. 2003) is ruled out to this confidence level; determining 
which individual assumptions are incorrect will required data beyond that of the current 
survey. These simulations (including a variety of other possible exoplanet distributions) are 
discussed in more detail in Nielsen et al. (in prep.) Nevertheless, our null detection in this 
survey sets strong upper limits on the distribution of young massive extrasolar planets >5 
AU from their primaries and provides valuable constraints for theories of planet formation 
and migration. 

4. Conclusions 

Wc obtained datasets for 54 stars (45 stars were observed in the southern sky at the 
VLT, 11 stars were observed in the northern sky at the MMT, and 2 stars were observed 
at both telescopes). In our VLT data, we achieved H band contrasts > 10 mag {5a) at a 
separation of 1.0" from the primary star on 45% of our targets and H band contrasts of > 
9 mag at a separation of 0.5" on 80% of our targets. With this degree of attenuation, we 
should be able to image (5(7 detection) a 5 Mjup planet 15 AU from the star around a 70 
Myr Kl star at 15 pc or a 5 Mjup planet at 2 AU from a 12 Myr M star at 10 pc. We believe 
that our SDl images are the highest contrast astronomical images ever made from ground 
or space for methane rich companions within 1" of their primary star. 

Eight tentative candidates were identified (none with S/N > 2 cr). Had these candidates 
been real, they would have possessed separations of 3 - 15.5 AU and masses of 2-10 Mj^p. 
However, none of the candidates were detected in second epoch observations. Thus, we find 
a null result from our survey. Nonetheless, our result still has serious implications for the 
distribution of extrasolar planets. In the course of our survey, we also discovered 5 new close 
stellar binary systems with measured separations of 0.14" to 0.26". 

For 20 of our survey stars, wc attained 50% completeness for 6-10 Mj„p planets at 
semi-major axes of 20-40 AU. Thus, our completeness levels are sufficient to significantly 
test theoretical planet distributions. From our survey null result, we can rule out (at the 
98%/2.0(T level) a model planet population using a constant distribution (N(a) oc constant) of 
planet semi-major axis out to a distance of 45 AU (a number of further exoplanet distribution 
models are considered in Nielsen et al. in prep). Our null detection in this survey sets strong 
upper limits on the distribution of young massive extrasolar planets >5 AU from their 
primaries and provides valuable contraints for theories of planet formation and migration. 

This pubhcation is based on observations made with the MMT and the ESO VLT at 
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Fig. 8. — 5a Contrasts for VLT SDI survey objects with H < 4.5 in the Fl (1.575 fim) 
filter. These contrast curves were generated by translating a 6x6 pixel (0.1"x0.r') box 
along a particular radial trajectory away from the center of the star and then calculating 
the standard deviation within that box as a function of radius. Curves were generated from 
the full reduced and differenced SDI data for each object (Fl(1.575 jjm) - F3a(1.625 /im) 
for two roll angles). 
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Fig. 9. — Same as Fig. [H] but for VLT SDI survey objects with 5.5 > H > 4.5. 
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Fig. 10. — Same as Fig. Obut for VLT SDI survey objects with 6.5 > H > 5.5. 
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Fig. 11. — Same as Fig. Obut for VLT SDI survey objects with 7.5 > H > 6.5. 
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Fig. 12. — Same as Fig. Obut for VLT SDI survey objects with H > 7.5. 
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Fig. 13. — 5(T contrasts for MMT SDI survey objects observed in May 2005. These contrast 
curves were generated by translating a 6x6 (0.1"x0.1") pixel box along a particular radial 
trajectory away from the center of the star and then calculating the standard deviation 
within that box as a function of radius. Curves were generated from the full reduced and 
differenced SDI data for each object (Fl(1.575 /im) - F3a(1.625 /im) for two roll angles). 
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Fig. 14. — Same as Fig. [13] but for MMT SDI survey objects observed in February 2006. 
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Fig. 15. — Seeing FWHM (averaged over each observation) vs. attained 5a contrast at 0.5" 
separation from the primary star for 10 of the stars presented in Fig. [9] with H magnitudes 
between 4.5 - 5.5. The error bars on seeing are the seeing variation (as measured by the 
standard deviation of the seeing) over each observation. For this sample of stars with roughly 
the same H magnitude, achievable contrast varies roughly inversely with the average seeing 
FWHM. Scatter in this plot is in part due to the fact that seeing FWHM can change 
considerably over a 20-40 minute long observation. 
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Fig. 16. — Minimum detectable planet separations for a 5 Mjup planet for the 10 objects 
in common between this survey and Masciadri et al. (2005) who used VLT NACO without 
SDI. For the purpose of comparison, we have adopted ages from Masciadri et al. (2005); 
we note our preferred age on the figure where our adopted ages differ from Masciadri et al. 
(2005). We translated theoretical 5 Mjup planet models (Burrows et al. 2003, Baraffe et 
al. 2003) into H magnitudes for these 10 cases then determined the minimum separation 
at which such a companion could be detected (at the 5a level) in our survey. For the 10 
objects in common between the surveys, our SDI survey attains lower minimum separations 
for 8 out of 10 objects (we note also that the two objects for which we did not attain lower 
separations were particularly low quahty AO/SDI datasets). 
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Fig. 17. — Planet detection completeness contour plots for a set of 4 typical program stars 
(upper left: AB Dor, upper right: DX Leo, lower left: GJ 182, lower right: GJ 799B). For a 
given mass and semi-major axis, 10,000 planets are simulated by our Monte Carlo method, 
over the expected distributions of eccentricity, orbital phase, and viewing angle. Given the 
parameters of the target star and the models of Burrows et al. (2003), we determine what 
fraction of the simulated planets are detectable at the 5a level given the contrast plot for that 
star. The contours show this detection probability across the 100,000 different combinations 
of mass and semi-major axis considered in this plot. The strong upper limit in mass is set by 
our conservative <1400 K limit for the methane break required for a robust SDI detection. 
In these models, we simply do not allow an object with Teff >1400 K to be detected, when 
in reality SDI can detect such non-methane objects (e.g. AB Dor C, Close et al. 2005b, 
Nielsen et al. 2005). For a complete set of planet detection completeness contour plots, see: 
http: / /exoplanet. as. arizona.edu/~lclose/SDI. html 
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Fig. 18. — Our 50% completeness levels. Combining the results of Fig. [TTl we consider 
individual values of the semi-major axis across the planetary mass range, and at each com- 
bination calculate the total number of stars in our survey (out of a total of 54) where the 
fraction of such planets, given by the Monte Carlo Simulation, that can be detected at the 
5a level is 50% or greater. Clearly, our survey is best able to place constraints on planets 
between 6 and 10 Mjup, and with semi-major axis between 20 and 40 AU. The decrease in 
sensitivity for masses >7 Mjup is due to the fact that such high mass planets are too hot 
to possess significant methane absorption if they are very young and, thus, are not ideal 
SDI targets. The higher completeness for 7-8 Mj^p planets for semi-major axis of 30 AU vs. 
semi-major axis of 40 AU is due to the small field of view of the SDI device; planets with 
semi-major axes > 30 AU can fall outside the SDI field in some of these cases. 
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Fig. 19. — A complete reduced dataset of AB Dor A (28 minutes of data at a series of 
rotator angles - 0°, 33°, 33°, 0°) from the VLT SDI device. Simulated planets have been 
added every 0.2" from the star (0.4", 0.6", 0.8", 1.0", 1.2", 1.4", 1.6", 1.8", 2.0", and 2.2") 
with AFl(1.575/im) = 9 mag (upper left, (attenuation in magnitudes in the 1.575 /im Fl 
filter), 10 mag (upper right), 11 mag (lower left) and 12 mag (lower right) fainter than the 
star. The 0.4" object falls within the inner dark circle (dark circle radius of 0.5", 0.5", 0.7", 
and 1.3" respectively for the 9, 10, 11, and 12 mag objects); the 2.2" object falls outside 
the frame aperture in a number of dither images and thus is detected with lower S/N than 
the other objects. These simulated planets are scaled from unsaturated images of AB Dor 
A taken right before the example dataset (and have fluxes and photon-noise in each filter 
appropriate for a T6 object). 
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Fig. 20. — Maximum achievable planet contrast {5a detection) vs. separation for 28 minutes 
of VLT SDI data for AB Dor A. To determine the maximum achievable planet contrast as 
a function of separation, we inserted and then attempted to retrieve simulated planets with 
a variety of separations and AFl contrasts appropriate for T5, T6, and T8 spectral types. 
The residual SDI noise curve for AB Dor A is also overplotted; the two curves agree well, 
giving us confidence in our measured contrast limits. 
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Fig. 21. — Maximum achievable H band planet contrast (Sa detection) vs. separation for 
28 minutes of VLT SDI data for AB Dor A. To determine the maximum achievable planet 
contrast as a function of separation, we inserted and then attempted to retrieve simulated 
planets with a variety of separations and AFl contrasts appropriate for T5, T6, and T8 
spectral types. AFl contrasts were converted to AH magnitudes using the magnitude offsets 
calculated in section 3.3. 
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Fig. 22. — Minimum detectable planet mass {5a detection) vs. separation (AU) for 28 
minutes of VLT SDI data for AB Dor A. To determine the minimum detectable planet mass 
as a function of separation, we inserted and then attempted to retrieve simulated planets with 
a variety of separations and AFl contrasts appropriate for T5, T6, and T8 spectral types. 
API contrasts were converted to AH magnitudes using the magnitude offsets calculated in 
section 3.3 and were then converted to absolute H magnitudes using the 2MASS apparent H 
magnitude and the Hipparcos distance for each star. Absolu te H magnitudes were converted 
into planet masses using the models of iBurrows et al.l (120031 ) and adopting a range of system 
ages from 50 - 100 Myr. For AB Dor, we should be able to image {5a detection) a 5 Mjup 
planet 12 AU from the star. For a complete set of minimum detectable planet mass vs. 
separation curves, see: |http://exoplanet.as.arizona. edu/^lclos e/SDI.html[ 
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Fig. 23. — Comparison of SDI contrast curve with other methods. The Lyot curve is for 
the 3.6m AEOS telescope (Hinkley et al. 2007) and the NICMOS curve coronograph curve 
is from HST (Schneider et ah 2003); otherwise curves are all from >8m class telescopes. 
We use the LQ Hya contrast curve from Masciadri et al. (2005) because this star (K2, 18 
pc vs. Kl, 15 pc) is the closest match from that work to AB Dor A (our SDI comparison 
star.) The SDI contrast curve has been converted from AFl=1.575/xm to AH contrasts 
appropriate for a T8 spectral type object. Inside 0.4", SDI contrasts are derived from the 
1-trajectory SDI contrast plot of AB Dor A; outside of 0.4", SDI contrasts are derived from 
our in-depth planet simulation case study of AB Dor A. For methanated companions, SDI 
provides improved contrast by 1-4 mag within 1" as compared to other methods. Past 1", 
narrowband imaging becomes less efficient and broad-band techniques (such as ADI; Marois 
et al. 2006) reach higher contrasts. 
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Fig. 24. — Minimum detectable mass vs. separation for a set of 4 typical program stars 
(upper left: AB Dor, upper right: DX Leo, lower left: GJ 182, lower right: GJ 799B). We 
convert our contrast curves in Amag units (from Figs.[8]tq 12|) in to minimum detectable mass 
vs. separation (in AU) using the models of lBurrows et al.l (l2003l ) and the distance to the star. 
To characterize the possible planets we expect to detect around each star, we simulated an 
ensemble of 10^ possible planets per star, assuming distributions for mass, eccentricity, and 
semi-major axis based on known radial velocity planets, as well as distributions for orbital 
phase and viewing angle. When combined with the properties of the individual target star 
and its measured contrast curve, we can determine what fraction of these simulated planets 
we expect to detect at the 5cr level (shown above each plot with the name of the target 
star). The ensemble of simulated planets is shown as small dots for each star in; simulated 
planets which are detected with the contrast attained by SDI are plotted in blue and those 
that remain undetected are plotted in red. Assuming each star possesses exactly one planet, 
we can assign a detection probability for that star from the percentage of simulated planets 
detected. For our 48 program stars which possess contrast curves, the average detection 
probability is 8.0%, the median detection probability is 4.1%, and the maximum detection 
probability is 47%. For GJ 799B (12 Myr M star at 10 pc), we can detect (at 5a) 5 Mj^p 
planet at 2 AU. 
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Fig. 25. — Expected number of planets detected. By taking the results of our Monte Carlo 
simulations, and assuming that each program star possesses exactly one planet, we can assign 
a detection probability for that star from the percentage of simulated planets detected. By 
adding these detection fractions for each star, we can compute the expected number of planets 
detected from our survey. We order the target stars by decreasing detection probability, and 
plot the total number of planets expected to be detected as a function of the number of stars. 
Over the entire survey, we expect to detect 3.8 planets, a 2(7 null result. Thus, our assumed 
distribution for the frequency (1 planet per star, hence 100%), semimajor axis distribution 
(N(a) oc constant), and luminosities (Burrows et al. 2003) of extrasolar planets is excluded 
at the 98% level by our extrasolar planet survey null result. 
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8.1 


5.3 


2.7% 




GJ 182 


04 


59 


34. 


,8 


+01 47 00.7 


26.7 


MlVe' 


0.016 


Li from FMS97 


10 


6.5 


15% 


very tentative planet candidate 
(4.S AU, -4 Mj^p) 
not detected at second epoch 


HIP 23309 


06 


00 


47. 


1 


-57 15 25.5 


26.3 


MO/1'" 


0.012 


Beta Pic 


10 


6.4 


16% 




AB Dor 


05 


28 


44, 


,8 


-65 26 54.9 


14.9 


KlIII" 


0.07 


AB Dor mg 


6.9 


4.8 


7.4% 


0.16" binary AB Dor C, 
Close et al. 2005 


GJ 207.1 


05 


33 


44. 


8 


+01 56 43.4 


16.8 


M2.5e" 


0.1 


Lowrance et al. 2005 


9.5 


7.1 


4.0% 




UY Pic 


05 


36 


56. 


,8 


-47 57 52.9 


23.9 


KOV° 


0.07 


AB Dor mg 


8 


5.9 


8.1% 




AO Men 


06 


18 


28. 


2 


-72 02 41.4 


38.5 


K6/7" 


0.012 


Beta Pic 


9.9 


7 


2.5% 


very tentative planet candidate 
(14 AU, ~4 Mj^p) 
not detected at second epoch 


HIP 30030 


06 


19 


08. 


1 


-03 26 20.0 


52.4 


GOVP 


0.03 


Tuc/Hor 


8 


6.6 


0.2% 




AB Pic 


06 


19 


12. 


,9 


-58 03 16.0 


45.5 


K2V'' 


0.03 


Tuc 


9.1 


7.1 


6.8% 


planetary mass companion, 
Chauvin et al. 2005 
very tentative planet candidate 
(15.5 AU, ~5 Mj„p) 
not detected at second epoch 


SRXl 


06 


22 


30. 


,9 


-60 13 07.1 


23.5 


GlV 


0.07 


AB Dor 


6.5 


5.2 


4.8% 




HD 48189A 


06 


38 


00. 


,4 


-61 32 00.2 


21.7 


G1/G2V'= 


0.07 


AB Dor 


6.2 


4.7 


1.8% 


0.14" binary 


BD +23 1978 


08 


36 


55. 


,8 


+23 14 48.0 


41.6 


K5V1 


0.035 


Montes ct al. 2001 


8.7 


6.5 






TT"! UMa 


08 


39 


11. 


,7 


+65 01 15.3 


14.3 


G1.5V'> 


0.21 


Li from WSH03 


5.6 


4.3 


0.1% 




LQ Hya 


09 


32 


25. 


,6 


-11 11 04.7 


18.3 


KOV" 


0.013 


Li from WSH03 


7.8 


5.6 


32% 




DX Leo 


09 


32 


43. 


,7 


+26 59 18.7 


17.7 


KOV 


0.115 


Her /Lyra 


7 


5.2 


1.8% 


very tentative planet candidate 
(2.6 AU, -10 Mj^p) 
not detected at second epoch 


TWA 22 


10 


17 


26. 


,9 


-53 54 28.0 


22 


M5^ 


0.01 




14 


8.1 






HD 92945 


10 


43 


28. 


,3 


-29 03 51.4 


21.6 


KIV" 


0.07 


AB Dor 


7.8 


5.8 


6.8% 


very tentative planet candidate 
(10.4 AU, ~6 Mj^p) 
not detected at second epoch 


GJ 417 


11 


12 


32. 


,4 


+35 48 50.7 


21.7 


GOV"' 


0.115 


Her /Lyra 


6.4 


5 


0.0% 




TWA 4 


11 


22 


05, 


,3 


-24 46 39.6 


46.7 


K4V* 


0.01 




9.1 


5.8 




0.78" binary 


TWA 25 


12 


15 


30, 


,8 


-39 48 42.0 


44.1 


MOS 


0.01 


TW Hydra 


11 


7.5 


14% 




RX J1224. 8-7503 


12 


24 


47, 


,3 


-75 03 09.4 


24.2 


K2'' 


0.016 


Li from AKCWM95 


11 


7.8 


13.5% 




RX J1231.9-7848 


12 


31 


56, 


,0 


-78 48 36.0 


50 


Ml" 


0.01 


Li from AKCWM95 


14 


9.6 






BK Dra 


14 


39 


00. 


2 


+64 17 30.0 


33.9 


GO" 


0.07 


AB Dor 


7.6 


6 


0.67% 


binary, 

Metchev & Hillenbrand 2004 


HD 135363 


15 


07 


56. 


,3 


+76 12 02.7 


29.4 


G5V'> 


0.0032 


Li from WSH03 


8.7 


6.3 


5.2% 


0.26" binary 


KW Lup 


15 


45 


47, 


,6 


-30 20 55.7 


40.9 


K2V^ 


0.002 


Li from NB98 


9.4 


6.6 


3.9% 




HD 155555AB 


17 


17 


25, 


5 


-66 57 04.0 


30 


G5IV+KOIV/V'" 


0.012 


Beta Pic 


7.2 


4.9 


0.06% 




HD 155555C 


17 


17 


27. 


,7 


-66 57 00.0 


30 


M4.5'" 


0.012 


Beta Pic 


13 


7.9 


47% 




HD 166435 


18 


09 


21. 


,4 


+29 57 06.2 


25.2 


GO" 


0.1 


RHK from Wright+ 04 


6.8 


5.4 


1.1% 





Table 1 — Continued 



T'.^ v.<^-f>4- 
larget 






Distance(pc) 


SpT 


Age(,GyrJ 


Age rtei 


V 


[jd 


Detect ability 




Comments 


HD 172B5BA 


18 45 26.9 


-64 52 16.5 


30 


A5IV/V'' 


0.012 


Beta Pic 


4.8 


4.3 


5.8% 






CD-64 1208 


18 45 37.0 


-64 51 44.6 


29.2 


K7" 


0.012 


Beta Pic 


10 


6.3 


9.9% 




0.18" binary 


HD 181321 


19 21 29.8 


-34 59 00.5 


20 


G1/G2V^ 


0.16 


Li from WSH03, 


6.5 


5 


0.09% 


very 


tentative planet candidate 














RHK from Gray+ 06 










(7 AU, ~5 Mj„p) 






















not 


detected at second epoch 


HD 186704 


19 45 57.3 


+04 14 54.6 


30.3 


GO'' 


0.2 


RHK from Wright+ 04 


7 


5.6 


0.0% 






GJ 799B 


20 41 51.1 


-32 26 09.0 


10.2 


M4.5e" 


0.012 


Beta Pic 


13 




36% 






GJ 799A 


20 41 51.2 


-32 26 06.6 


10.2 


M4.5e" 


0.012 


Beta Pic 


11 


5.2 


18% 






GJ 803 


20 45 09.5 


-31 20 27.1 


9.94 


MOVe° 


0.012 


Beta Pic 


8.8 


4.8 


33% 


very 


tentative planet candidate 
























(3 AU, ~2 Mj„p) 






















not 


detected at second epoch 


HIP 112312A 


22 44 57.8 


-33 15 01.0 


23.6 


M4e8 


0.012 


Beta Pic 


12 


7.2 


34% 


very 


tentative planet candidate 
























(6.2 AU, -8 Mj^p) 






















not 


detected at second epoch 


HD 224228 


23 56 10.7 


-39 03 08.4 


22.1 


K3V 


0.07 


AB Dor 


8.2 


6 


5.8% 






More Distant Young Stars 


TWA 14 


11 13 26.5 


-45 23 43.0 


66.7 


MO" 


0.01 


TW Hydra 


13 


8.7 


7.8% 






RXJ 1243.6-7834 


12 43 36.7 


-78 34 07.8 


150 


MO" 


0.008 


Li from AKCWM95 


13 


8.7 






0.068" binary 


Stars with known RV planets 


eEri 


03 32 55.8 


-09 27 29.7 


3.22 


K2V' 


0.8 


Benedict et al. 2006 




1.9 


0.1% 




Kellner et al. 2007, 
























Janson et al. 2007 


HD 81040 


09 23 47.1 


+20 21 52.0 


32.6 


GOV 


2.5 


Li from SUZT06 


7.7 


6.3 


0.0% 




1 


HD 128311 


14 36 00.6 


+09 44 47.5 


16.6 


K0<1 


0.63 


RHK from Gray+ 03 


7.5 


5.3 


0.0% 




1 


Nearby Solar Analogues 


HD 114613 


13 12 03.2 


-37 48 10.9 


20.5 


G3V" 


4.2 


Li from RGP93, 


4.8 


3.3 


0.0% 


















RHK from Gray+06 












HD 201091 


21 06 53.9 


+38 44 57.9 


3.48 


K5Ve" 


2.0^ 


RHK from Gray+ 06 


5.2 


2.5 


0.0% 






eind A 


22 03 21.7 


-56 47 09.5 


3.63 


K4.5V" 


1.3 


Lachaume et al. 1999 


4.7 


2.3 


0.09% 




Geifiler et al. 2007 


GJ 862 


22 29 15.2 


-30 01 06.4 


15.4 


K5Vo" 


6.3^ 


RHK from Gray+ 06 


7.7 


5.3 


0.0% 







"'"In general, wc have only determined Ca R'HK ages for stars with spectral types Kl or earlier, but in the case of these two K5 stars, 
we have only the R'HK measurement on which to rely for age determination. The calibration of Mt. Wilson 

S-index to R'HK for K5 stars {B-V '^1.1 mag) has not been well-defined {Noyes et al. 1984; specifically the photospheric subtraction), and hence 

applying a R'HK vs. age relation for K5 stars is unlikely to yield useful ages. Although we adopt specific values for the ages of these stars, it would be more accurate to state simply that these 
stars have ages 

>1 Gyr. As a result, almost all simulated planets are too faint to detect around these stars, so the precise error in the age does not significantly affect our final results, 
^derived from the Hipparcos survey. Ferryman et al.(1997) 

'^ages for stars with cluster memberships from Zuckerman and Song (2[)[)4), otherwise, ages are either lithium ages, calcium R,HK ages, or an average of both. 
Acronyms for lithium and calcium age references: AKCWM95: Alcala, Krautter, Schmitt, Covino, Wichmann, and Mundt 1995, 
FMS97: Favata, Micela, and Sciortino 1997, LH98: Li and Hu 1998, NB98: Neuhauser and Brandner 1998, 
RGP93: Randich, Gratton, and Pallavicini 1993, SUZT06: Sozetti, Udry, Zucker, Torres, Beuzit, et al. 2006, 
WSH03: Wichmann, Schmitt, and Hubrig 2003 

'^from the CDS Simbad service 

'^from the 2MASS Survey, Cutri et al. (2003) 

* Spectral reference key: e: Houk and Cowley 1975, f: Houk and Cowley 1988, 
g: Zuckerman and Song 2004, h: Li and Hu 1998, i: Cowley, Hiltncr, and Witt 1967. 

j: Christian and Mathioudakis 2002, k: Leaton and Pagel 1960, 1: Favata, Barbera, Micela, and Sciortino 1995, 



m: Zuckerman, Song, Bessell, and Webb 2001, n: Gliesse and Jahreifi 1991, o: Houk 1978, 
p: Cutispoto, Pallavicini, Kuerster, and Rodono 1995, 

q: Montes, Lopez- Santiago, Galvez, Fernandez-Figueroa, De Castro, and Cornide 2001, 

r: Bidelman 1951, s: Zuckerman, Webb, Schwartz, and Becklin 2001, t: Houk and Smith-Moore 1988, 

u: Alcala, Krautter, Schmitt, Covino, Wichmann, and Mundt 1995, v: Houk 1982, 

w: Gliesse and Jahreifi 1979, x: Henry Draper Catalog 
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Table 2. VLT SDI Observation Log 



yj Ui KjK^ t 




DTT 


NDTT 


TVi'i'f^l T-f^iYT* ( TYllTm'i'PQ 1 
XVJtCti I ilililU.tC& J 


HIP 1481 


2004-11-15 


14 


6 


56 




2005-11-24 


16 


5 


26.7 




2005-11-25 


16 


5 


26.7 




2005-11-27 


16 


5 


26.7 


ERX6 


2004-11-14 


22 


4 


29.3 




2004-11-16 


22 


4 


58.7 


ERX8 


2004-11-17 


22 


4 


58.7 


HIP9141 


2004-09-27 


14 


6 


56 


BD +05 378 


2005-02-01 


32 


3 


25.6 


HD 17925 


2003-08-14 


7.5 


16 


40 




2003-08-16 


4 


30 


40 




2003-08-17 


4 


30 


20 




2004-02-02 


1 


120 


20 




2004-11-16 


4.1 


17 


46.5 




2004-11-17 


4.1 


17 


46.5 


LH98 062 


2004-02-03 


14 


9 


21 


eEri 


2004-09-19 


0.6 


160 


64 


V834 Tau 


2005-01-25 


10 


9 


24 




2005-02-01 


10 


9 


24 


GJ 182 


2004-02-02 


7 


17 


39.7 




2005-11-22 


20 


4 


26.7 




2005-11-24 


20 


4 


26.7 




2005-11-27 


20 


4 


26.7 


HIP23309 


2005-01-30 


24 


4 


25.6 




2005-01-31 


24 


4 


51.2 


AB Dor 


2004-02-02 


5 


24 


20 




2004-09-28 


12 


7 


28 




2004-11-16 


10.4 


8 


27.7 


GJ 207.1 


2005-01-27 


32 


3 


25.6 


UY Pic 


2004-11-16 


14 


6 


28 




2004-11-17 


14 


6 


56 
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Table 2 — Continued 



Object Date DIT NDIT Total Exp (minutes) 



ivien 


onn/L no it? 

ZUU^t-UZ-UO 


1 4 


Q 


91 
Zl 




onnc; 11 11^ 

ZUUO-ii-iO 


QO 


i 


17 c; 




9nn'i-1 1-94 


30 


1 

X 


10 

xu 


Ad r IC 


ofin/i 11 1/1 
ZUU4-ii-i4 


9n 
zU 


/I 
4 


9fi 7 
ZO. 1 




90(14 11 1 ^ 

ZUU^~ 1 1- 1 (J 


90 

ZU 




Of, 7 

ZU. ( 




zuuo- i i-ZZ 


90 
ZU 


/I 
4 


1 Q Q 

io.o 




9nnc; 1 1 9^; 
zuuo- i i-ZO 


90 
ZU 


/I 
4 


Oo.o 




9nn4 11 18 


1 9 
xz 


7 


98 
zo 




9nn4-i 1-1 Q 


1 9 


7 


98 




9004-1 1-1 7 


fi 5 


1 1 

± X 


93 8 




9004 11 18 


U.U 


1 1 

X X 


zo.o 


RD -lO-^ 1078 


9005-01-97 


94 


4 


9^ fi 




2005-01-28 


24 


4 


25 6 


LO Hva 


2004-02-02 


5 


24 


40 




2004-1 9-08 


1 4 


VJ 


28 




9004- 1 9-1 4 

^yjyj t: X Z( ± t: 


1 4 


U 


98 


DX Leo 


2004-02-05 


3 


38 


19 




2005-12-04 


14 


6 


28 




2005-12-19 


14 


6 


28 


i WA ZZ 


onnc; ni oc; 
zUUo-Ui-zo 


QO 
OZ 


i 


48. 


HD92945 


2004-02-05 


5 


24 


60 


TWA 14 


2005-01-28 


32 


3 


25.6 




2005-01-29 


32 


3 


25.6 


TWA 4 


2004-02-02 


7 


17 


9.92 


TWA25 


2005-01-28 


32 


3 


25.6 


RX J1224.8-7503 


2004-02-02 


40 


3 


20 




2005-01-16 


30 


3 


60 




2005-01-27 


30 


3 


120 


RX J1231.9-7848 


2004-02-05 


20 


6 


20 


RXJ 1243.6-7834 


2004-02-02 


5 


24 


40 


HD 114613 


2004-02-02 


1 


120 


40 
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Table 2 — Continued 



Object Date DIT NDIT Total Exp (minutes) 



KW Lup 


2004-09-15 


22 


4 


14.7 




2004-09-16 


24 


4 


22.7 




2004-09-17 


24 


4 


24 


HD155555AB 


2003-08-14 


7.5 


16 


10 




2003-08-15 


7.5 


16 


20 




2003-08-16 


7.5 


16 


10 




2003-08-17 


7.5 


16 


10 




2004-09-16 


10 


9 


30 




2004-09-18 


14 


6 


28 


HD155555C 


2003-08-14 


30 


4 


40 




2003-08-16 


30 


4 


40 


HD172555A 


2003-08-17 


5 


24 


20 




2004-09-17 


5 


15 


25 




2004-09-18 


5 


15 


6.25 




2004-09-19 


5 


15 


18.8 


CD-64 1208 


2003-08-17 


20 


6 


40 




2004-09-16 


15 


6 


30 


HD181321 


2003-08-15 


7.5 


16 


40 




2004-09-18 


11 


8 


29.3 


GJ799B 


2003-08-16 


20 


6 


40 




2003-08-17 


20 


6 


30 




2004-09-19 


15 


6 


30 


GJ799A 


2003-08-16 


20 


6 


40 




2004-09-16 


10 


9 


30 




2004-09-19 


15 


6 


30 


GJ803 


2003-08-14 


7.5 


18 


56.2 




2003-08-15 


10 


12 


40 




2003-08-17 


7.5 


16 


40 




2004-09-17 


6 


15 


30 




2001-09-18 


10 


9 


30 



eIndA 2004-09-18 0.5 192 



48 
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Table 2 — Continued 



Object Date DIT NDIT Total Exp (minutes) 



GJ 862 


2003-08-15 


10 


12 


40 




2003-08-16 


10 


12 


40 




2004-09-19 


13 


7 


48.2 


HIP112312A 


2001-09-19 


25 


1 


66.7 


HD224228 


2003-08-16 


10 


12 


40 




2003-08-17 


20 


6 


40 




2004-10-08 


14 


6 


28 




2004-10-20 


21 


4 


28 
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Table 3. MMT SDI Observation Log 



Object 


Date 


DIT 


NDIT 


Total Exp (minutes) 


V577 PerA 


2006-02-12 


20 


7 


37.3 




2006-02-13 


21.5 


7 


40.1 


HIP30030 


2006-02-12 


30 


5 


30 


TTi UMa 


2006-02-13 


5.8 


13 


40.2 


HD 81040 


2006-02-12 


11.7 


13 


40.3 


LQ Hya 


2006-02-12 


8 


19 


40.5 


DX Leo 


2005-05-01 


10 


13 


34.7 


GJ 417 


2005-04-30 


7 


17 


31.7 


HD 128311 


2006-02-12 


4 


19 


60.8 


EK Dra 


2005-05-01 


20 


7 


37.3 


HD 135363 


2005-05-01 


30 


5 


40 


HD 166435 


2005-04-30 


7 


17 


31.7 




2005-05-01 


'/' 


17 


31.7 


HD 186704 


2005-05-01 


10 


13 


17.3 


HD 201091 


2005-04-30 


20 


7 


37.33 
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Table 4. Limiting H mag {5a) at 0.5" 



Object 


API 


Separation(AU) 


AH (T8 SpT) 






e Eri 


y.4 ± U.iz 


i.oi 


1 n o 
iU.O 


1 o o 

iz.z 


1/17 

i4. t 


e incl A 


iU.O ± U.iz 


i.oi 


ii.O 


1 o o 
io.o 


1 p. 
io 


rlJJ zUiUyi 


o.Uo ± U.oz 


1 7/1 

i. /4 


o.yo 


1 1 

ii.O 


1 Q Q 

io.o 


TTT~V 1 1 /I CI Q 

rlJJ ii4Dio 


D.io ± U.ZD 


1 n o 
lU.z 


1 


1 n 
iU.o 


Q 7/1 

o. /4 


rlU i / yzo 


y.by ± U.14 


o.iy 


iU.D 


1 /I o 

i4.5 


1/17 

i4. / 


rlJJi ( zoooA 


n 1 /I 1 n 1 o 

y. 14 ± u.iz 


ID 


1 n 
iU 


1 /I Q 

i4.o 


1 1 n 
ii.y 


TTi UMa 


O.U4 ± U.iO 


/ .i4 


Q m 

o.yi 


1 Q o 

io.z 


T O /I 

iz.4 


XJT\ /I o 1 on A 


o.04 ± U.UOZ 


T n Q 
iU.o 


y.4i 


T /I O 

i4.z 


iz.O 


Tono 
(orJoUo 


n /I _i_ n nm 

y.o4 ± u.uyi 


/I n'T 
4.y 1 


1 n /I 
iU.4 


1 o 
io.z 


1 o 
io.z 


ArS JJor 


n n/i _i_ n ni n 

y.U4 ± u.uiy 


'7 A'7 


n m 

y.yi 


1 /I Q 

i4.o 


ion 

io.y 


riJJioooooAr) 


o.o/ ± U.i4 


10 


R 7/1 

D. /4 


1 1 

ii.D 


o oi 
y.zi 




/. li) ± U.zo 


inn 

iu.y 


O.DD 


1 Q 7 

io. / 


1 o 
iz 


rlJJioiozi 


/.4z ± U.lo 


1 n 
iU 


Q on 
o.zy 


1 Q Q 

io.o 


1 1 Q 
ii.O 


ortAi 


y.yo ± u.u/y 


117 
ii. / 


1 n Q 
iU.o 


1 R 

ib 


1/11 

i4.i 


Lrj / yyA 


( .4o ± U.Uoz 


1 1 

O.ii 


O.OO 


io.o 


io.o 


UA Leo 


o.z4 ± u.iy 


Q Q7 
O.O / 


nil 

y.ii 


1/1/1 
i4.4 


1 Q O 

io.z 




y.oi ± U.ZO 


7 70 
/ . /Z 


1 n /I 
iU.4 


1 f; 7 
io. / 


^ A Q 
i4.o 


Voo4 iau 


y.Uo ± U.io 


^; 7/1 
0. /4 


n nF; 

y.yt) 


1 f; Q 
io.o 


^ A p. 
i4.0 


riU iOD4oo 


o.4z ± U. i / 


iz.O 


n on 

y.zy 


1/17 

i4. / 


1 O 7 

iz. I 


LLi -tiya 


n CO 1 nip 
y.oz ± U.iD 


n 1 7 

y.i / 


1 n 7 
iU. / 


1 Q 
io.o 


1 c: 
io 


rlJJ ioD/U4 


T 1 _L n nm 

I.L6 ± u.uyi 


1 1 
io.i 


Q 

o 


1 Q 
io.D 


1 1 o 
ii.z 


HD92945 


9.91 ± 0.0099 


10.8 


10.8 


16.6 


14.9 


UY Pic 


9.96 ± 0.11 


11.9 


10.8 


16.7 


14.8 


HD224228 


9 ± 0.15 


11 


9.87 


15.9 


14.2 


EK Dra 


7.85 ± 0.39 


17 


8.72 


14.7 


12 


HIP 1481 


9.22 ± 0.13 


20.5 


10.1 


16.3 


13.2 


CD-64 1208 


9.33 ± 0.087 


14.6 


10.2 


16.5 


14.2 


HD 135363 


7.9 ± 0.27 


14.7 


8.77 


15.1 


12.8 


HIP23309 


8.45 ± 0.092 


13.1 


9.32 


15.7 


13.6 


GJ 182 


8.01 ± 0.16 


13.3 


8.88 


15.3 


13.2 


V577 PerA 


8.9 ± 0.33 


16.9 


9.77 


16.2 


13.6 
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Table 4 — Continued 



Object 


AFl 


Separation (AU) 


AH (T8 SpT) 




Mh 


HIP9141 


8.92 ± 0.29 


21.2 


9.79 


16.3 


13.2 


HIP30030 


6.91 ± 0.17 


26.2 


7.78 


14.4 


10.8 


KW Lup 


8.76 ± 0.091 


20.5 


9.63 


16.3 


13.2 


ERX8 


9.4 ± 0.2 


18.6 


10.3 


17.2 


14.4 


ERX6 


9.38 ± 0.4 


24.6 


10.2 


17.1 


13.6 


AO Men 


6.91 ± 0.33 


19.2 


7.78 


14.8 


11.9 


AB Pic 


9.65 ± 0.027 


22.8 


10.5 


17.6 


14.3 


GJ 207.1 


7.5 ± 0.094 


8.41 


8.37 


15.5 


14.4 


HIP112312A 


9.09 ± 0.27 


11.8 


9.96 


17.1 


15.2 


BD +05 378 


8.31 ± 0.088 


20.3 


9.18 


16.4 


13.4 


TWA25 


9.5 ± 0.035 


22 


10.4 


17.9 


14.7 


RX J1224.8-7503 


7.16 ± 0.024 


12.1 


8.03 


15.9 


14 


HD155555C 


10.5 ± 0.085 


15 


11.4 


19.3 


16.9 


TWA 14 


8.38 ± 0.03 


33.3 


9.25 


18 


13.9 
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Table 5. Limiting H mag {5a) at 1.0" 



Object 


AFl 


Separation(AU) 


AH (T8 SpT) 






e Eri 


LL.O ± O.Z 


00 
o.zz 


100 
Iz.z 


1/11 
i4.i 


io.o 


e Ind A 


1 _i_ n 1 c 
iz ± U.iD 


0.00 


ion 

iz.y 


1 
iO.z 


17/1 

i lA 


rlJJ zUiUyi 


y.4z ± U.Uo 




1 n 
iU.o 


1 Q 

iz.o 


1 1 
io.i 


TJT\ 1 1 /I CI Q 

rlJJ ii4Dio 


/.z4 ± U.io 


on K 
zU.o 


Q 1 1 
O.ii 


ii.o 


n nA 

y.y4 


rlU i / yzo 


ii.o ± o.iy 


1 n y( 

iU.4 


100 

iz.z 


iD.4 


ICO 

io.o 


TTT^ 1 '70KKK A 

rlJJi ( zoooA 


ii.Z ± u.uyo 


Qn 
oU 


101 
iz. i 


1 /I 

iO.4 


1 A 
i4 


TTi UMa 


y.zo ± U.i4 


1 /I Q 

i4.c) 


1 n 1 
iU.i 


1/1/1 
i4.4 


1 Q C 

io.o 


UT\ /I Q1 Qfl A 


y.o/ ± U.z4 


01 ^7 

zl. i 


1 n 7 
iU. / 


1 /I 
iO.4 


1 Q 7 

io. 1 


LrJoUo 


^ Ci ^ _i_ n no 
iO. 1 ± O.Oo 


n n /I 

y.y4 


11c 

ii.o 


1 /I 
iO.4 


1 C /I 

iO.4 


ArS JJor 


ii ± U.LI 


1 /I n 

i4.y 


1 1 n 

ii.y 


1^7 
io. i 


1 Q 
io.o 


riJJioooooAr) 


/.o ± U.U4D 


Qn 
oU 


Q 1 7 

o.i / 


io.i 


1 n 7 
iU. / 




5.44 ± U.Uo 


zi. / 


n 1 

y.oi 


^ A 
i4.o 


10c 
iz.D 


rlJJioiozi 


o.Do ± U.U45 


on 
zU 


n K 

y.o 


i4.D 


io.i 


DrtAi 


ii.Z ± U.io 


OQ K 

Z6.0 


101 

iz.i 


I 7 Q 

II .0 


1 K A 

iO.4 


T'7nn A 

Lrj ^ yyA 


y.OO ± U.i4 


1 n 
iU.z 


iU.4 


io.o 


1 c 

io.o 


UA Leo 


y.yo ± u.uoy 


17 7 


1 n Q 
iU.o 


1 a 

io 


1 /I Q 

i4.o 


T c(;o 

LjtJ OOZ 


1 n V _i_ n 10 
iU. / ± U. iz 


1 /I 
iO.4 


1 1 

ii.o 


1 n 

io.y 


1 

io 


Voo4 iau 


1 n _i_ n 10 
iU.Z ± U.io 


1 Q 
id.O 


111 

i i. i 


1/^/1 

iO.4 


1 p; 7 
io. 1 


JJT\ 1 /^/^ /I Q 

i004o0 


n no _i_ n n^i 
y.yo ± U.Uoi 


ot^ 

zo.z 


1 n Q 
iU.o 


1^0 
iO.z 


1 /I 

i4.z 


T r\ Wtto 


ii ± U.Uoo 


ICQ 

io.o 


1 1 n 
ii.y 


17 c: 
i / .0 


1 
iD.z 


nU ioD/U4 


/.OO ± U.Uoz 


on Q 
oU.o 


Q 00 

o.zz 


io.o 


11/1 
ii.4 


HD92945 


10.8 ± 0.062 


21.6 


11.7 


17.5 


15.8 


UY Pic 


11.5 ± 0.033 


23.9 


12.4 


18.3 


16.4 


HD224228 


10.8 ± 0.11 


22.1 


11.7 


17.7 


16 


EK Dra 


8.86 ± 0.14 


33.9 


9.73 


15.7 


13 


HIP 1481 


10.8 ± 0.046 


41 


11.7 


17.9 


14.8 


CD-64 1208 


9.88 ± 0.54 


29.2 


10.8 


17.1 


14.8 


HD 135363 


8.65 ± 0.025 


29.4 


9.52 


15.8 


13.5 


HIP23309 


10 ± 0.051 


26.3 


10.9 


17.3 


15.2 


GJ 182 


10.2 ± 0.15 


26.7 


11.1 


17.6 


15.5 


V577 PerA 


10 ± 0.062 


33.8 


10.9 


17.4 


14.8 
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Table 5 — Continued 



Object AFl Separation(AU) AH (T8 SpT) m^^ M 



HIP9141 


10.5 ± 0.028 


42.4 


11.4 


1^ 




14.9 


HIP30030 


8.3 ± 0.09 


52.4 


9.17 


15, 


.8 


12.2 


KW Lup 


9.86 ± 0.17 


40.9 


10.7 


17, 


.3 


14.2 


ERX8 


10.7 ± 0.12 


37.1 


11.6 


18, 


.5 


15.7 


ERX6 


10.6 ± 0.12 


49.3 


11.5 


18, 


.4 


14.9 


AO Men 


7.9 ± 0.015 


38.5 


8.77 


15, 


.8 


12.9 


AB Pic 


10.8 ± 0.013 


45.5 


11.7 


18, 


.8 


15.5 


GJ 207.1 


8.74 ± 0.089 


16.8 


9.61 


16, 


.8 


15.7 


HIP112312A 


10.6 ± 0.068 


23.6 


11.5 


18, 


.7 


16.8 


BD +05 378 


9.52 ± 0.074 


40.5 


10.4 


17, 


.6 


14.6 


TWA25 


10.5 ± 0.18 


44.1 


11.4 


18, 


.9 


15.7 


RX J1224.8-7503 


8.04 ± 0.16 


24.2 


8.91 


16, 


.8 


14.9 


HD155555C 


10.8 ± 0.043 


30 


11.7 


19, 


.6 


17.2 


TWA 14 


8.74 ± 0.047 


66.7 


9.61 


18, 


.3 


14.2 



Table 6. Star/Planet Projected Minimum Detectable Separations for 5 and 10 Mj^p Planets 



Ubject 


Age (Myrj 


Distance (pc) 


Separation 5 Mj^p 


(AU) 


beparation iU Mjup 


/ A TT^ 

(AU) 








Burrows et al. 


Baraffc ct al. 


Burrows et al. 


Baraffc ct 


AB Dor 


70 


14.94 


13.45 


20.92 


6.28 


10.31 


AO Men 


12 


38.48 


30.01 


* 


11.54 


20.01 


BD+05 378 


12 


40.54 


16.21 


28.37 


7.70 


9.32 


CD -64 1208 


12 


34.21 


11.29 


17.45 


5.82 


8.55 


DX Leo 


115 


17.75 






16.68 


33.54 


EK Dra 


70 


33.94 










eEri 


800 


3.22 










eind A 


1300 


3.63 








3.81 


GJ 174 


160 


13.49 






12.01 


17.54 


GJ 182 


12 


26.67 


6.67 


16.27 


4.80 


5.60 


GJ 207.1 


100 


16.82 






10.77 


15.48 


GJ 417 


115 


21.72 










GJ 799A 


12 


10.22 


1.94 


4.60 


1.12 


1.64 


m 7QQR 


1 9 






1 QA 


119 


1 19 


GJ 803 


12 


9.94 


1.39 


2.39 


1.09 


1.09 


GJ 862 


6300 


15.45 










HD 114613 


4200 


20.48 










HD 128311 


630 


16.57 










HD 135363 


3 


29.44 


8.54 


11.78 


3.24 


6.77 


HD 155555 AB 


12 


30.03 






25.23 




HD 155555 C 


12 


30.03 


3.30 


3.30 


3.30 


3.30 


HD 166435 


100 


25.24 






29.03 





al. 



Table 6 — Continued 



Object 


Age (Myr) 


Distance (pc) 


Separation 5 Mj„p 


(AU) 


Separation 10 Mj„p 


(AU) 


nU i 1 ZDDO A 


iz 


on OQ 


on 1 T 
zU.l r 


OS 

Z5.D0 
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14. oz 


XJT^ 1 vnott 
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zUUU 












TTr\ OO/IOOQ 

rlJJ ZZ'iZZo 


/ U 


oo nc 
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y.ou 
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/U 


zo.oU 


on 7Q 




1 n 
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1 no 
io.yo 


TTT^ /I SI SO A 


7n 
lU 


01 fiV 
zi.D / 






1 n OQ 

ly.yo 


o7 nc; 
z/ .yo 


rlJJ oiU4U 


zoUU 


Qo 










rlU oOOo 


Qn 
oU 


/in on 
4y.zy 


OT 1 1 
Z ( .11 


fi7 C^O 
D ^ .OZ 


1 1 Q/l 

11. o4 


oo 1 s 
zz. lo 


rlU ioO(U4 


zUU 


Qn ofi 










rlJJ yUo4 


on 


Q'T 1 C; 


1 n QO 
ly.dz 


Q1 on 


V n^ 


n on 

y.zy 


rlJJ yzy4o 


/ U 


01 c;7 
zl.o / 


O/i 

z4.oy 


A A nn 
44. UU 


R n/i 
D.U4 


1 n 1 /I 
iU.i4 


WTP 1 1 0Q1 A 
jnir iiZoiZ A 


1 

±z 


OQ R1 
Zo. Di 


O.04 


1 Q 

o.iy 


Z.DU 


BQ 


TTTTD 1 /I Q1 

rlir i4oi 


Qn 
oU 


/in nc; 
4U.yo 


Qn 71 
oU. / 1 


KCi QQ 

oy.oo 


1 o (^n 
iz.Dy 


on QQ 
zU.oo 


TTTT3 oQonn 


iz 


OR o<5 
ZD. ZD 


R QQ 


1 no 
iD.Uz 


o on 

z.sy 


A nn 

4.yy 


HIP 30030 


30 


52.36 






53.93 




AB Pic 


30 


45.52 


19.12 


40.05 


9.56 


13.65 


HIP 9141 


30 


42.35 


37.70 


78.78 


15.25 


22.45 


KW Lup 


2 


40.92 


6.14 


7.36 


4.50 


4.50 


LQ Hya 


13 


18.34 


3.30 


4.95 


2.02 


2.75 


RXJ1224.8-7503 


16 


24.17 


6.53 


17.16 


3.63 


4.59 


TWA 14 


10 


66.67 


14.00 


27.33 


7.33 


10.67 


TWA 25 


10 


44.05 


9.25 


14.10 


4.85 


7.93 



Table 6 — Continued 



Object 


Age (Myr) 


Distance (pc) 


Separation 5 Mj„p 


(AU) 


Separation 10 Mj„p 


(AU) 


UY Pic 


70 


23.87 


19.81 


29.12 


5.73 


11.46 


V577 Per A 


70 


33.77 


81.73 




17.90 


27.02 


TTi UMa 


210 


14.27 











* — means that such an object is too low in mass to be detected with our current survey contrast 
level for that star 



-63- 



Table 7. Binary Properties 



Object Separation Position Angle 



SDI survey discoveries 



AB Dor AC^ 


0.16" 


127° 


HIP 9141 


0.15" 


355° 


HD 48189AC 


0.14" 


143° 


HD 135363 


0.26" 


132° 


CD -64 1208 


0.18" 


95° 


SDI survey confirmations 


RXJ 1243.6-7834^ 


0.068" 


1717351° 


LH 98 062 


2.4" 


354° 


TWA 4 


0.78" 


3° 


EK Dra 


0.67" 


176° 



^Separation and position angle from Close et al. 2005b. 
For updated photometry and astrometry see Close et al. 
2007b. 

^As RXJ 1243.6-7834 is nearly an equal-magnitude bi- 
nary, we were unable to determine which star was the pri- 
mary (as selected by Brandncr et al. (2000)) and thus 
present two values for the position angle (assuming each 
star is the primary in turn). 
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